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In light-front dynamics, the regularization of amplitudes by traditional cutoffs imposed on the 
transverse and longitudinal components of particle momenta corresponds to restricting the integra- 
tion volume by a non-rotationally invariant domain. The result depends not only on the size of 
this domain (i.e., on the cutoff values) , but also on its orientation determined by the position of the 
light-front plane. Explicitly covariant formulation of light front dynamics allows us to parameterize 
the latter dependence in a very transparent form. If we decompose the regularized amplitude in 
terms of independent invariant amplitudes, extra (non-physical) terms should appear, with spin 
structures which explicitly depend on the orientation of the light front plane. The number of form 
factors, i.e., the coefficients of this decomposition, therefore also increases. The spin-1/2 fermion 
self-energy is determined by three scalar functions, instead of the two standard ones, while for the 
elastic electromagnetic vertex the number of form factors increases from two to five. In the present 
paper we calculate perturbatively all these form factors in the Yukawa model. Then we compare the 
results obtained in the two following ways: (i) by using the light front dynamics graph technique 
rules directly; (ii) by integrating the corresponding Feynman amplitudes in terms of the light front 
variables. For each of these methods, we use two types of regularization: the transverse and longitu- 
dinal cutoffs, and the Pauli-Villars regularization. In the latter case, the dependence of amplitudes 
on the light front plane orientation vanishes completely provided enough Pauli-Villars subtractions 
are made. 

PACS numbers: ll.10.-z, 11.25.Db, 13.40.Gp 
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I. INTRODUCTION 

Light-Front Dynamics (LFD) is extensively and suc- 
cessfully applied to hadron phenomenology, relativis- 
tic few-body systems, and field theory. For reviews 
of theoretical developments and applications see, e.g., 
Refs. [HQ. In the LFD framework, non-perturbative ap- 
proaches to field theory were developed in Refs. [1, 0, Q 
and in Refs. [1, 0> Q • I n spite of some essential differ- 
ences, these approaches proceed from the same starting 
point, namely, they approximate the state vector of the 
system by a truncated one. The problem is then solved 
without any decomposition in powers of the coupling 
constant. Doing that, one should carry out the renor- 
malization procedure non-perturbatively. This is a non- 
trivial problem which is at the heart of ongoing intense 
research. However, before renormalization, one should 
regularize the amplitudes, both in the perturbative and 
non-perturbative frameworks. The explicit dependence 
of these amplitudes on the cutoffs is not unique. It is 
determined by the method of regularization. 

It is of utmost importance to understand the origin and 
the implications of this dependence if one wants to ad- 
dress the question of non-perturbative renormalization. 
As we shall show in the present article, this is the only 
way to identify the structure of the counterterms needed 
to recover full rotationally invariant renormalized ampli- 
tudes in LFD. The question of the non-perturbative de- 
termination of the counterterms in truncated Fock space 



is discussed in Ref. [8J. 

The rules for calculating amplitudes can be derived di- 
rectly, either by transforming the standard T-ordering 
of the S-matrix into the ordering along the Light Front 
(LF) time (the LFD graph technique rules or from 
quantized field theory on the LF plane 0. An alterna- 
tive method consists in expressing well-defined Feynman 
amplitudes through the LF variables and integrating over 
the minus-components of particle momenta [{J, [13] ■ Such 
an approach was applied to the derivation and study 
of the LF electromagnetic amplitudes [lj, [12] and to 
the analysis of different contributions to the electromag- 
netic current, resulting from the LF reduction of the 
Bethe-Salpeter formalism A study of electroweak 

transitions of the spin-1 mesons, based on using the 
LF plane of general orientation was carried out in 
Refs. [H], [Hi- The comparison of perturbative ampli- 
tudes obtained from the LFD graph technique rules with 
those derived from the Feynman approach is, in general, 
not trivial, as shown in Ref. [161 ]. 

Concerning the regularization of amplitudes in LFD, at 
least two important features should be mentioned. First, 
if a given physical process is described by a set of LF di- 
agrams, each partial LF amplitude usually diverges more 
strongly than the Feynman amplitude of the same pro- 
cess. The statement holds true regardless of the origin 
of the LF contributions: either from the rules of LFD, 
or from the Feynman amplitude. This increases the sen- 
sitivity of the result to the choice of the regularization 
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procedure and may be a source of the so-called treacher- 
ous points [l7|. Second, the LF variables (and, hence, the 
integration domain with the cutoffs imposed on it) explic- 
itly depend on the LF plane orientation, which means the 
loss of rotational invariance. Because of this extra depen- 
dence, standard decompositions of such regularized LF 
amplitudes into invariant amplitudes are not valid. The 
total number of invariant amplitudes (and the number 
of form factors which are the coefficients in this decom- 
position) increases, as compared to the case when the 
rotational invariance is preserved. For example, the LF 
electromagnetic vertex (EMV) of a spin-1/2 particle is 
determined by five form factors rather than by two. In 
order to cancel the extra contributions (which depend on 
the LF plane orientation) one needs to introduce in the 
interaction Hamiltonian new specific counterterms. 

This complication is especially dramatic in non- 
perturbative approaches, where the Fock space trunca- 
tion is another source of the rotational symmetry vio- 
lation. A given Feynman diagram may generate a few 
time-ordered ones with intermediate states containing 
different number of particles. When they are truncated, 
the rotational invariance is lost even for invariant cut- 
offs. The interlacing of the two sources of the violation 
of rotational symmetry makes non-perturbative analysis 
of the counterterm structure extremely involved. One 
should therefore separate to a maximal extent the prob- 
lems coming from the regularization procedure and from 
the Fock space truncation. That can be effectively done 
in the explicitly covariant formulation of LFD within the 
perturbative framework in a given order in the coupling 
constant. 

In the present paper we study in detail this problem 
for the spin-1/2 fermion perturbative self-energy and the 
elastic EMV in the Yukawa model within the framework 
of explicitly covariant LFD [l[ . The latter deals with the 
LF plane of general orientation uj-x = u>ot — uj-x = 0, 
where w is a four-vector with uj 2 = 0. In the particular 
case uj = (1,0,0,-1) we recover the standard approach 
on the plane t + z = 0. Due to w, which is transformed as 
a four-vector under rotations and Lorentz boosts, we can 
keep manifest rotational invariance throughout the calcu- 
lations. Dependence of amplitudes on the LF plane orien- 
tation turns now in their dependence on the four-vector 
to. The latter participates in the construction of the spin 
structures in which the regularized initial LF amplitude 
can be decomposed, on equal footing with the particle 
four-momenta. This generates extra (^-dependent) spin 
structures with corresponding scalar coefficients (e.g., 
electromagnetic form factors). The number of the extra 
spin structures and their explicit forms are determined by 
general physical principles (more precisely, by the parti- 
cle spins and the symmetries of the interaction), i.e., they 
are universal for any model and do not depend on par- 
ticular features of dynamics. Whereas, the dependence 
of the extra coefficients on particle four-momenta is de- 
termined by the model. This allows to separate general 
properties, related to LFD itself, from model-dependent 



effects. 

We shall proceed in the following two ways, both for 
the self-energy and the EMV. In the first way, we calcu- 
late these quantities by the LFD graph technique rules, 
taking into account all necessary diagrams. In the sec- 
ond way, we start from the standard Feynman amplitudes 
and integrate them in terms of the LF variables. For both 
ways, we have to introduce cutoffs on the LF variables, 
and this fact already implies the contribution of extra 
(w-dependent) structures and their corresponding form 
factors. The regularized self-energy and the EMV cal- 
culated by means of the LFD graph technique rules do 
not coincide in general with their counterparts obtained 
from the Feynman amplitudes. For the EMV case, the 
vertex found from the Feynman amplitude with the cut- 
offs imposed on the LF variables, also differs from that 
calculated in a standard way, by the Wick rotation with a 
spherically symmetric cutoff or by the Pauli-Villars (PV) 
regularization. All these differences disappear when we 
deal with integrals which are finite from the very begin- 
ning (e.g., due to the PV regularization). They disappear 
also in the renormalized amplitudes, though renormal- 
ization procedures (counterterms, etc.) are drastically 
different for different regularization schemes. 

Within covariant LFD, the perturbative QED self- 
energy and the EMV in the channel of the fermion- 
antifermion pair creation have been studied earlier [l8| . 
The main subject was to extract the physical (u>- 
independent) contributions from the corresponding am- 
plitudes and to renormalize this physical part only. The 
present analysis is devoted to a more detailed treatment 
of the self-energy and the EMV. We calculate both phys- 
ical and non-physical contributions in the two ways men- 
tioned above and investigate the influence of the regular- 
ization procedure on the whole amplitudes and on their 
subsequent renormalization. The Yukawa model which 
we use reflects some features of QED but it is simpler 
from the technical point of view. 

The paper is organized as follows. In Sec. [Til we briefly 
describe the LFD graph technique rules and apply them 
to calculate the fermion self-energy. We use two different 
regularization procedures, the transverse and longitudi- 
nal LF cutoffs or PV subtractions, either for the bosonic, 
or simultaneously for both bosonic and fermionic propa- 
gators. Then we calculate the fermion self-energy, start- 
ing from the manifestly invariant Feynman amplitude 
expressed through the LF variables and regularized in 
the same way as the LFD one. We compare the results 
obtained in both approaches and analyze how they are 
affected by the choice of regularization. In Sec. Hill we 
repeat analogous steps for the fermion EMV. For this 
purpose, we derive the LF interaction Hamiltonian which 
includes fcrmion-boson and fcrmion-photon interactions. 
We then construct the complete set of the LF diagrams 
which contribute to the EMV. We use again the non- 
invariant LF cutoffs and the invariant PV regularization. 
The LFD form factors are compared to those obtained in 
terms of the Feynman amplitude with the same type of 
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regularization. General discussion of our results is pre- 
sented in Sec. HVl Sec. [V] contains concluding remarks. 
The technical details of some derivations are given in the 
Appendices. 



sions for the corresponding amplitudes read 



9 



(2tt) s 



6{wk)5(k 2 - f/)<Tk 



x {4 + m)6{u-q)8{q A - m J )(fq 



II. 



THE FERMION SELF-ENERGY 



The fermion self-energy is the simplest example of how 
an extra spin structure is generated by rotationally non- 
invariant cutoffs in LFD. To make the situation more 
transparent, we will calculate the self-energy indepen- 
dently in the two following ways: (1) by applying the 
covariant LFD graph technique rules; (2) by using the 
four-dimensional Feynman approach. In each case we 
consider two different types of regularization of diver- 
gent integrals: either the traditional rotationally non- 
invariant cutoffs or the invariant PV regularization. We 
then renormalize the amplitudes and compare the results 
obtained within these two methods. 



A. Calculation in light-front dynamics 

1. Light-front diagrams and their amplitudes 

We calculate in this section the fermion self-energy in 
the second order of perturbation theory, using t he g raph 
technique rules of explicitly covariant LFD [l|, 0, [t9( . We 
do it in details in order to explain the rules on a concrete 
example. 
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FIG. 1: Two contributions to the LFD fermion self-energy 
— E(p): the two-body intermediate state (a) and the contact 
term (b). The solid, dashed, and dotted lines represent, re- 
spectively, the fermion, the boson, and the spurion. Here 
p = pi — ior\ , where lot\ is the four-momentum attached to 
the initial (or final) spurion line. See text for the explanation. 



The self-energy is determined by the sum of the 
two diagrams shown in Fig. [U 



E(p) = S 2 b(p) + S /c (p). 



(1) 



They correspond to the two-body contribution and the 
fermion contact term, respectively. Analytical expres- 



s /c(p) 



q) -, (2a) 



7- - iO ' 

2 ,,2\ j4i 



(2tt) 3 J 2uj{p-k) 



6{uj-k)5{k 2 - ^ 2 )d 4 k 



(2b) 



where g is the coupling constant of the fermion-boson in- 
teraction, m and [i are the fermion and boson masses, re- 
spectively. In the covariant LFD graph technique, all the 
four-momenta are on the corresponding mass shells. This 
is due to the fact that the propagators are proportional to 
delta- functions: 8(ui-k)S(k 2 —iJ 2 ) is the boson propagator, 
(4+ m)8(ui-q)5(q 2 — m?) is the fermion one. Each theta- 



function 9{lo-1) selects only one value of Iq = y 1 + mf , 
of the two possible ones allowed by the corresponding 
delta-function 6(l 2 — mf). There is no any conservation 
law for the components of particle four-momenta in the 
direction of u> (or for the minus-components, in the stan- 
dard version of LFD). The conservation is restored by 
the spurion four-momentum lot which enters the delta- 
function (p + lot — k — q) . The factor l/(r — iO) is the 
spurion propagator and the r-integration is performed 
in infinite limits. To avoid misunderstanding, we em- 
phasize that spurions are not true particles and do not 
affect particle counting. They serve as a convenient way 
to describe the departure of intermediate particles off the 
energy shell. The term "spurion" itself is used for shorter 
wording only. For this reason, the intermediate state in 
the self-energy (|2a[) contains one fermion and one scalar 
boson only. The self-energy is supposed to be off-energy- 
shell, i.e., n ^ 0. 

Integrating by means of the delta- functions over d^q, 
dr, and dfco, we get 

g 2 f {$- ybr + m)6[u)-(p - k)\ d 3 k 



S/ C (P) 



(2tt) e 



2lo-(p — k)-i 



d 3 k 



l± 

(27t) 3 7 2w(p-*)2ek' 



where £k = ko 



and 



i 2 - (p - kf 
2lo-{ P - k) 



2e k ' 
(3a) 

(3b) 



(4) 



Let us go over to the LF variables. First, we denote 
x = (uj-k)/ (co-p) (equivalent to k + /p + in standard non- 
covariant LFD on the surface t + z = 0). We then split 
the three- vector k into two parts: k = k^ + kn , which 
are, respectively, perpendicular and parallel to the three- 
vector ui. Since S(p) is an analytic function of p 2 , we 
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may calculate it for p 2 > 0, while its values for p 2 < p = 0. Using the kinematical relations from AppendixlAl 
are obtained by the analytical continuation. If p 2 > 0, we one can rewrite Eqs. © as 
may perform our calculation in the reference frame where 



^2b(p) 
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16tt 3 



d 2 k\ 



m)dx 



32tt 3 (wp) 



d 2 k x 



+ m 2 x - 

+°° dx 



p 2 x(l-x)+n 2 (l-x) 32tt 3 (w-p) 



d 2 k. 



dx 



x(l — x) ' 



x(l — x) 



Both T,2b(p) and S/ C (p) are expressed through integrals 
which diverge logarithmically in x and quadratically in 
|kj_|. Possible regularization procedures are discussed 
below. 

Following Ref. 0], we will use the matrix representa- 
tion 



^2b{p) 



A(p 2 ) + B(p 2 )1 + C(p 2 )^ 
m uj-p 



9 2 C f , 



mi/) 
uo-p ' 



(6a) 
(6b) 



where the coefficients A, B, and C are scalar functions 
which depend on p 2 only. They are independent of uj. 
The coefficient Cf c is a constant. The self-energy is thus 
obtained by summing up Eqs. (J6j> : 



E(p) = g 2 A( P 2 ) + B(p 2 )^- + [C( P 2 ) + C fc 



mt/i 
uj-p 



(J) 



(5a) 
(5b) 



r 



Note that in the expression ([7]) an additional spin struc- 
ture proportional to i/j appears, as compared to the stan- 
dard four-dimensional Feynman approach. 



2. Regularization with rotationally non-invariant cutoffs 



In order to regularize the integrals over d 2 k± , we intro- 
duce a cutoff Aj_ , so that k\ < . Since some integrals 
over dx diverge logarithmically at x = and /or x — 1 , we 
also introduce (where it is needed) an infinitesimal posi- 
tive cutoff e, assuming that x may belong to the intervals 
e < x < 1 — e and 1 + e < x < +co. The correspond- 
ing analytical expressions for the functions A(p 2 ), B(p 2 ), 
C(p 2 ), and Cf c were found in Ref Q: 



A( P 2 ) 
B( P 2 ) 

C(P 2 ) 
C fc 



16tt 2 J 



m 



dk i / dx 



16tH J 
1 



f 



dk i / dx 



k]_ + m 2 x — p 2 x(l — x) + /i 2 (l — x) 
l-x 



l-e 



k^ + m 2 x — p 2 x(l — x) + /i 2 (l — x) ' 



dki 



dx 



k\ + m 2 — p 2 (l — x 



(1 — -^[k^ + m 2 x — p 2 x(l — x) + /i 2 (f — x) 



32ir 2 m 



dkj 



dx 



dx 



x(l — x) Ji+ e x(l — x) 



(8a) 
(8b) 
(8c) 
(8d) 



We imply in the following that A^ 3> max{|p 2 |, m 2 , /i 2 } 
and e <C 1. The dependence of physical results on the 
cutoffs is eliminated by taking the limits A_l — > oo, e — > 0. 



Retaining in Eqs. ([5]) all terms which do not vanish in 
these limits, we get 
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o, to , A± to f 1 \m 2 x - p 2 x(\ - x) + ^ 2 (1 - x)~\ ,„ . 

„, «. m , A to f 1 , . . , \m 2 x — p 2 x(l — x) + u 2 (l — x)l , , , 

S P ) = - ttt-j log — + -TpTl, / log ^ 1 , (9b) 

l07T z TO 107T Z Jq [ TO 

,9, 1 to 2 — u 2 A i I/99 9 m\ , , 

C ^ 3 ) = ~ 32^to" ^ 7 - m£ - 32^ (- 2 - * - 2 ^ ^ -J , (9c) 



We have not integrated over dx in Eqs. (|9a|l and (|9bj) 
because the results of the integrations are rather long. It 
is interesting to note that C(p 2 ) does not depend on p 2 . 
Due to this, we will denote in the following C(p 2 ) = C — 
const and, for shortness, C = C + Cf c . 

Each of the two quantities, C and C/ c diverges like 
A 2 L log(l/e). The strongest divergencies cancel in the 
sum C = C + Cf c , but the latter differs from zero and, 
moreover, has no finite limit when Aj_ — * 00: 



self-energy is the w-dependence of the integration do- 
main. A standard way free from this demerit is the use of 
the PV regularization, since in that case the cutoffs have 
no more relation to u>. In the language of LFD, the PV 
regularization consists in changing the propagators as 



6(wk)S(k 2 - v 2 ) -> d(uj-k)[5(k 2 - u 2 ) - S(k 2 



■A) 



~ to 2 - /i 2 Aj_ 

C = -TT- 2 log 

l07T z TO TO 



1 



32tt 2 to 



to 2 - 11 2 - 2u 2 log ■ 



M 
(10) 

Since the two diagrams shown in Fig. [1] exhaust the full 
set of the second-order diagrams which contribute to the 
fcrmion self-energy, we might expect the disappearance 
of every u>- dependent contribution from the amplitude. 
As far as this does not take place, the only source of the 
cj-dependence is the use of the rotationally non-invariant 
cutoffs for the LF variables and x. Indeed, if we write 
these variables in the explicitly covariant form 



for scalar bosons, and 



2 _ (^-fc)(M 2 

K _L — L P 

UJ-p 



uj-k 
ijj-p 



OJ-k 

LO-p' 



it becomes evident that both of them depend on lo. 
Introducing the cutoffs A^ and e, we restrict an in- 
dependent integration domain, which inevitably brings 
w-dependence into the regularized quantities. We will 
demonstrate this feature in more detail in Sec. HVI bv us- 
ing a very simple and transparent example. Note that 
without adding new, aj-dependent, counterterms in the 
interaction Hamiltonian, this dependence is not killed 
by the standard renormalization. The renormalization 
recipe must be therefore modified 0- 



3. Invariant Pauli-Villars regularization 

As we learned above, the source of the appearance of 
the extra (w-dependent) term in the regularized fermion 



-I- m)6(uj-q)S(q 2 — to 2 ) — > 

(wg) [{4+m)5{q 2 - to 2 ) - {4+ mi )8{q 2 - to 2 )] 



for fermions. This procedure is equivalent to introducing 
additional particles (one PV fermion with the mass m\ 
and one PV boson with the mass n±), whose wave func- 
tions have negative norms. If needed, more subtractions 
can be done till all integrals become convergent. After 
the calculation of the integrals and the renormalization, 
the limits [i\ — > 00 and toi — > 00 should be taken. In the 
case of the fermion self-energy, we may regularize either 
the boson propagator only, or the fermion one, or both 
simultaneously. Hereafter we will supply PV-regularized 
quantities with the superscript " PV, b" (when only the 
bosonic propagator is modified) or " PV, b + f" (when 
both bosonic and fermionic propagators are modified). 

Let us now calculate the PV-regularized coefficients 
A(p 2 ), B{p 2 ), and C. We can start from the expres- 
sions ([9]), in spite of their dependence on the "old" cutoffs 
Aj_ and e. Indeed, the integrals in Eqs. ([5]) become reg- 
ular, provided enough PV subtractions have been made. 
If so, they have definite limits at Aj_ — > 00 and e — > 0. 

The integrals for A(p 2 ) and B(p 2 ) become convergent 
after the regularization by a PV boson only: 



J 
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A pv > b (p 2 ) = A(p 2 ,m >f j,)-A(p 2 ,m,n 1 ) = 
B pv ' b (p 2 ) = B{p 2 ,m,fi)-B{p 2 ,m,^) = 



m 
16tt 2 

m 
16tt 2 



<ix log 



m 2 x — p 2 x(l — x) + /i 2 (l — x) 
m 2 x — p 2 x(\ — x) + /if (1 — a;) 



dx (1 — x) log 



p 2 x(l — x) + /t 2 (l - x) 
p 2 x(l — x) + ^1(1 - x) 



r 



(lla) 
(lib) 



The situation differs drastically for the coefficient C. 
After the bosonic PV regularization we get 



C 



PV,b 



1 



32tt 2 to 



(P 2 



Mr) 



1 + 2 log ■ 



A i 



2/i 2 log- -2 M 2 log- 



Mi 



(12) 



Since the result is still divergent for Aj_ — > oo, this reg- 
ularization is not enough. The additional fcrmionic PV 
regularization requires some care because C is a coeffi- 
cient at the spin structure mi/)/ (uj-p) which itself depends 
on m. Hence, one should regularize the quantity mC: 



mC 



(mCX 



= mC{m, /i) — mC(m, /ii) 
miG{m\, /i) + toiC(toi, /Ui) = 0. 



(13) 



We see that after the double PV regularization the extra 
structure in Eq. l[7]). proportional to i/, disappears, as it 
should. Note that Eq. (fl~3|) holds for arbitrary (i.e., not 
necessary infinite) PV masses m\ and 



In order to calculate A re n(p 2 ) and B ren (p 2 ) we can use 
the initial functions A(p 2 ) and B(p 2 ) regularized either 
by the non- invariant cutoffs, Eqs. (|9a|) and (|9b[). or by 
means of the PV regularization, Eqs. (fTTj) . Any choice 
leads to the same result: 



Aren (,P ) 



#ren(p 2 ) 



where 

4>i{x) 



m 
16tt 2 

m 
16tt 2 



dx [</>i(x) - </> 2 (a;)] , 



(18a) 



dx (1 - x) [01 (x) + <fo(aO] 



(18b) 



m 2 x - p 2 x(l - x) + /i 2 (l - x) 



^ m 2 x 2 + /i 2 (l — x) 

2m 2 x(2-3x + x 2 ) 
m 2 x 2 + /i 2 (l — x) 



The remaining integrations over dx in Eqs. (|18j) are sim- 
ple but lengthy. 



^. Renormalization procedure 

The renormalizcd self-energy is obtained by using the 
standard procedure: 



^ b+ f{p) = ^ PV ' b +f(p)-c 1 -c 2 ( I (-m), (14) 



where 



Cl 



u(p)E py ' b +/(p)w(p) 



2m 



1 f gg^W(p) 



Finally, 



Ke V n b+f ( P )=9 2 



A (p 2 ) + B ren (p 2 )- 



(15a) 
(15b) 



(16) 



with 

A ren (p 2 ) = A(p 2 ) - A(m 2 ) + 2m 2 [A'{m 2 ) + B'(m 2 )}, 



(17a) 



B ren {p ) = B{p z ) - B{m 2 ) - 2m z [A\m 2 ) + B'(m )] 



(17b) 



B. Calculation in the four-dimensional 
Feynman approach 

1. Regularization with rotationally non-invariant cutoffs 

We showed above that the regularized fermion self- 
energy calculated within LFD with the traditional trans- 
verse (A j.) and longitudinal (e) cutoffs contains an extra 
spin structure depending on oj. In order to understand 
the reasons of this behavior, we calculate here the self- 
energy in another way, following Ref. Q. We start from 
the standard four-dimensional Feynman expression 



Ef(p) 



(2*)* 



d 4 k 



[k 2 - /i 2 + i0] [0 - k) 2 - m 2 + iO] ' 
' . (19) 

but perform the integrations in terms of the LF variables, 
with the corresponding cutoffs. For this purpose, we in- 
troduce the minus-, plus-, and transverse components of 
the four-momentum k: 



k- = ko - k z , k + = k + k„, k ± = (k x , k 



and analogously for p. As in Sec. Ill A 11 we take, for 
convenience, the reference frame where p = 0. In this 
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frame p_ = p 2 /p+. Denoting k + = xp + , we get 



Sjr(p) 



- — £ / d 2 k i / da; 
32tt 4 



' — oo 



[k-p+x — k 2 L —Li 2 + iO] 
— ft + to) 



[(p_ - fc_)p+(l - x) - - m 2 + iO] 



(20) 



with 



1 x 
¥ = 7p+ k - + 2^- p+ ~ 7 -L" k -L- 

The integral over dk- is calculated by using the principal 
value prescription: 



/+oo pL 
&_(...) = lim / dk_(. 



(21) 



This integral is well defined unless x = or x = 1 . If a; = 
or x = 1, the integral (|21j) diverges on the upper and 
lower limits. So, the infinitesimal integration domains 
near the points x — and x = 1 (the so-called zero 
modes) require special consideration. We introduce a 
cutoff e in the variable x and represent T,(p) as a sum 



S F (p) = £ p+a (p) + £ zm (p), 



(22) 



where S p+a (p) incorporates the contributions from the 
regions of integration over dx (— oo < x < oo), excluding 
the singular points x — and x — 1, while the zero- mode 
part E zm (p) involves the integrations in the e- vicinities 
of these two points. 

The calculation is carried out in Appendix [BJ After 
closing the integration contour by an arc of a circle (see 
Appendix |SJ , the integral for E p+a (p) is represented as 
a sum of the pole and arc contributions and has the form 



16tt 3 



d 2 k A 



(p 1 — ft + m)dx 



+ m 2 x — p 2 x(l — x) + ^i 2 (l — x) 



d 2 k t 



-t-oo 



dx 



327r 3 (w-p) J y i+e x(l — x) 

Comparing Eq. ([23]) with Eqs. we see that 

S p+a(p) = ^2b(p) + E/ C (p), 



(23) 



(24) 



where both terms on the right-hand side are regularized 
at x = and x = 1. The pole plus arc contributions to 
Eq. (|20[) reproduce the result given by the sum of the two 
LFD diagrams shown in Fig. [T] Hence, 



A(p 2 ) + B(p 2 )l + C p+a ^ 

m uj-p 



(25) 



where A(p 2 ) and B(p 2 ) are defined by Eqs. (|9a|) and (|9b|) . 
while C p+a coincides with C, Eq. (fTUjl . 

As mentioned above, the zero mode contribution re- 
sults from the divergence (at L — > oo) of the integral 
over fc_, which occurs when x = 0,1. This divergence is 
determined by the leading fc_-term in the numerator of 
Eq. (|2H)) . that is by ^7 + fc_. In explicitly covariant LFD, 
the matrix 7+ turns into The zero- mode contribution 
is therefore 



E zm (p) = g 2 C z 



mi/j 
oj-p 



(26) 



where C zm is calculated in Appendix [Bl and has the form 



C 2 



We thus find 



i« 2 lo s — 



(27) 



+ 



1 



327r 2 m 



/i 2 - 2^ 2 log ■ 



p+a 



0. 



(28) 



since C p + a — C is given by the right-hand side of 
Eq. |10]). Substituting Eqs. flS3) and (26]) into Eq. (|22|) . 
and taking into account Eq. {28]) , we finally get 



Sjr(p) = .g 2 



^b 2 ) 



m 



(29) 



We see that after the incorporation of the pole, arc, 
and zero- mode contributions, the w-dependent term in 
the self-energy disappears, even before the renormaliza- 
tion. Note that although we used the same cutoffs, the 
formula ([29]) does not coincide with the expression ob- 
tained by using the LFD graph technique rules, since the 
sum C(p 2 ) + Cf c is not zero [it is given by Eq. (fTU]) ]. One 
might think that the dependence of the self-energy on 
w is an artefact of the LFD rules, whereas the indepen- 
dence of the Feynman approach on to is natural, since we 
started with the Feynman expression (|19[) which "knows 
nothing" about u>. It is not so, since the initially di- 
vergent integral for £(p) acquires some sense only after 
regularization, and the latter has been done in terms of 
the cutoffs imposed on the LF variables. The indepen- 
dence of E_f(p) on u> looks thereby as a coincidence. We 
shall see in Sec. MI 51 that, in the EMV case, in contrast 
to the self-energy one, the LF cutoffs applied to the ini- 
tial Feynman integral do result in some dependence of 
the EMV on the LF plane orientation. 

After the renormalization, Eq. (|2"9"]l reproduces the ex- 
pression (|16p obtained earlier for the self-energy found 
within LFD and regularized by the invariant PV method. 



2. Invariant regularization 

We briefly recall in this section familiar results known 
from the standard four-dimensional Feynman formalism 
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which is completely independent from LFD. It may serve 
as an additional test of the results obtained above. 

The fermion self-energy is given by Eq. (|19p . It is con- 
venient to use the following decomposition: 



A F (p 2 )+B F {p 2 )^- 
m 



(30) 



where Af(p 2 ) and Bf(p 2 ) are scalar functions. Note 
that they do not coincide with A(p 2 ) and B(p 2 ) from 
Eq. (|29p . since we use here another regularization pro- 
cedure. Applying the Feynman parametrization, we can 
rewrite Eq. (fl"9)) as 



S F (p) = 



dx ■ 



19 

16tH 



d 4 k' 



a-x) 



ft' + m 



[k 1 



m 2 x + p 2 x(l — x) — ji 2 (l — x) + iO] 2 ' 



(31) 



where we introduced k' = k — xp. Going over to Eu- 
clidean space by means of the Wick rotation k' Q = ik' 4 
with real k\ and regularizing the divergent integrals by 
an invariant cutoff \k' | = k' + k'i < A 2 (assuming 



that A 2 > 



{m 2 ,n 2 , 



|p 2 |}), we can perform the four- 



dimensional integration and get 



Sf(p) 



x lo 



16tt 2 



dx [(1 — x)p" + m] 



A 2 



l x — p 2 x{\ — x) + fi 2 (l — x) 



- 1 



(32) 



The terms of order 1 /A 2 and higher are omitted. Com- 
paring the right-hand sides of Eqs. (f3"0"|) and (j52"|) , we find 



Af( P 2 ) = 
B F {p 2 ) = 



m 
16?r 2 

m 
32^ 



log — - 1 
m z 

, A 2 
log — - 1 



16tt 2 



dx log 



dx (1 — x) log 



m 2 x -p 2 x(l - x) + /i 2 (l - x) 



m 

m 2 x — p 2 x(l — x) + ^i 2 (l — x) 



(33a) 
(33b) 



r 



The renormalized self-energy Hp. ren (p) is found 
from Eq. (H}, changing E py ^+^(p) by S F (p). The 
corresponding renormalized functions ^_F, r en(p 2 ) and 



B F ,ren(p 2 ) coincide with A r 



1 

(p 2 ) and B ren (p 2 ) given 



by Eqs. (|18|) . We thus reproduce again Eq. <jl6|> for 

Sf, ren(p) ■ 

It is easy to verify that using the PV subtraction for 
the boson propagator (instead of the cutoff A 2 ), 

^ F V ' b (p) = ^F(p,m,fi) - £f(p,to,Mi)> 

leads, in the limit \i\ — > 00, to the same expression 
for the renormalized self-energy. We have therefore 

^PV,b 



tip) 



l (p)=E^(p). (34) 



To conclude, any of the considered ways of regulariza- 
tion (either with the non-invariant cutoffs A^ and e or 
with the single bosonic PV subtraction) can be used in 
order to get the correct expression for the renormalized 
self-energy in the Feynman approach, while the double 
(bosonic + fermionic) PV subtraction is required in the 
case of LFD. 



III. 



FERMION ELECTROMAGNETIC VERTEX 



We can now proceed to the calculation of the spin- 1/2 



fermion elastic EMV, 



T p , which is connected with the 



matrix element of the electromagnetic current J p by the 
relation 



J p = eu(p')T p u(p), 



(35) 



where e is the physical electromagnetic coupling con- 
stant, p and p' are the initial and final on-mass-shell 
fermion four-momenta (p' 2 — p 2 — m 2 ). Assuming P-, 
C-, and T-parity conservation, T p is defined by two form 
factors depending on the momentum transfer squared 
Q 2 = -{p' -p) 2 : 



U 1 pU = u 



Fi(Q 2 h P 



*F 2 (Q 2 ) 
2m 



-o pv q v 



(36) 



We omit for simplicity the bispinor arguments and de- 
noted q = p' -p, a pv = i{jpj v - j v j p )/2. 

We consider here, as in Sec. [Ill the Yukawa model 
which takes into account interaction of fermions with 
scalar bosons, while the EMV "dressing" due to fermion- 
photon interactions is neglected. The fermion-boson in- 
teraction is treated perturbatively, up to terms of order 
9 2 - 

At that order, the EMV must be renormalized. The 
standard renormalization recipe consists in the subtrac- 
tion r r p en = T p — Zi^f p with the constant Z\ found from 
the requirement u(p)T r p en u(p) = u(p)"f p u(p). This leads 
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to the following well-known expressions for the renormal- 
ized form factors: 

Fr n (Q 2 ) = 1 + Fi(Q 2 ) - Fi(0), F^ n (Q 2 ) = F 2 (Q 2 ). 

(37) 

We shall follow the same ideology that we exposed 
above for the self-energy: independent calculations of 
the EMV are performed, within covariant LFD and the 
Feynman approach, both for non-invariant and invariant 
rcgularization. However, in contrast to the self-energy 
case, the use of rotationally non-invariant cutoffs results 
in the appearance of new structures (and form factors) 
in the EMV, even if one starts from the standard four- 
dimensional Feynman expression. By this reason, the 
renormalization of the two physical form factors only, as 
prescribed by Eqs. (|37|) . is not enough to get the full 
rcnormalized EMV. 



that standing between the exponents acts on all the func- 
tions to the right of it. In momentum representation, the 
action of the operator 1 / (ito-d) on a function f(x) reduces 
to the multiplication of its Fourier transform f(k) by the 
factor l/(u>-k). 

It is easy to modify the Hamiltonian (f3"5)) in or- 
der to incorporate interactions between fermions and 
scalar bosons. The equation of motion for the Heisen- 
berg fermion field operator 'J, in the absence of scalar 
bosons, looks like (ift — m)^> = — ejp^ 1 where A is 
the Heisenberg photon field operator. If we introduce 
a scalar boson field the equation of motion becomes 
(ift — m)^ = (— eA\ — g$>)ty. Since the latter equation of 
motion is obtained from the previous one by the substi- 
tution A\ — > A\ + (g/e)&, it is enough to make the same 
substitution in the Hamiltonian (|38p . everywhere except 
in the operator 1/(iljj-Da) 0j- The Hamiltonian thus 
becomes 



A. Calculation in light-front dynamics 

1. Light-front interaction Hamiltonian involving 
electromagnetic interaction 

Before going over to the consideration of the EMV in 
covariant LFD, one should derive the interaction Hamil- 
tonian which involves both fermion-boson and fermion- 
photon interactions. Although all the diagrams calcu- 
lated below are generated by the graph technique rules 
formulated in Ref . [l| , this is another way to explain their 
origin. 

We will not give here a detailed derivation. The cor- 
responding procedure is exposed in Ref. Q. We derived 
there the LF Hamiltonian describing a system of interact- 
ing fermion and massless vector boson fields. This expres- 
sion holds also for a fermion-photon system. When the 
photon field is taken in the Feynman gauge, the Hamil- 
tonian in Schrodinger representation has the form 



H mt (x) = -e^4V + e 2 ^4 



2i(uj-D A ) 



4V, 



(38) 



where ip and A are the free fermion and photon fields, 
respectively, and the operator 1/(iu;-Da) in coordinate 
space acts on the coordinate along the four-vector to (for 
convenience we denote it for a moment as X- , since u) in 
standard LFD has only the minus-component): 



1 



f(x-) = cxp 



j-d 



-ujA 



~~7i i ex P 



iui-d 



ujA 



/(*-)} (39) 



and l/{iu>-d) is the free reversal derivative operator: 

j^gfi*-) =-\j d V-< x - - V-)f(V-)> (4°) 

where e(x) is the sign function. The operators l/(iu)-d) 
inside the exponents act on the functions luA only, while 



H mt {x)= - ij[gtp + e4]ip 



2iuj-Da 



[g<p + eJf\ ip. (41) 



We see that the LFD interaction Hamiltonian (j4Tj) in- 
volves also, besides the usual term — ip [gip + eJf\ ip de- 
scribing ordinary fermion-boson and fermion-photon in- 
teractions, the so-called contact terms which are nonlin- 
ear in the coupling constants. Note that if we expand 
the Hamiltonian (|4ip in powers of e, this expansion con- 
tains an infinite number of terms. Such a peculiarity is 
connected with the photon spin and with the gauge we 
have chosen. 

In this paper, we are not interested in studying elec- 
tromagnetic effects, but focus on the interaction between 
fermions and scalar bosons. We therefore restrict the 
Hamiltonian to the first order in the electromagnetic cou- 
pling constant e, neglecting the terms of order e 2 and 
higher. The result is 



H v 



where 



* =H 

H fbl 
H fb2 

-^eml 
H e m2 



fbl 



H 



m 



H„ r - 



H, 



-giptp(p, 
g 2 ipip 



2iu-d 

eipAilj, 



(pip, 



H, 



egtp ip 



1 



em 3 



2iuj-d 



A+A 



2iuj-d ' 



(42) 

(43a) 
(43b) 
(43c) 

<p)1>, (43d) 



2 \ iuj-o 



ito-d 



)A 



ito-d 



)A 



ioj-d 



(43e) 



The contact terms are Hfb2i Hem2, and H em3 . The op- 
erators l/(iwd) inside the squared brackets act on ujA 
only. 
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2. Light-front diagrams and their amplitudes 




FIG. 2: Triangle LF diagram 



Since the amplitude of the process which we are inter- 
ested in is proportional to eg 2 , we should collect together 
the matrix elements from the Hamiltonian (|42[) in the 
first, second and third orders of perturbation theory. It 
can be written schematically as 

< Hf bl H em i > + < HfblH em 2 > + < H em 3 > . 

Note that although the matrix element of the sec- 




UT 1 



FIG. 3: Pair creation by a photon. The double solid line 
stands for an antifermion 




p k P ' p k 



(a) (6) 
FIG. 4: Left (a) and right (b) contact terms 



non-covariant version of LFD). It comes from the analy- 
sis of the pure scalar "EMV" (i.e., when all the particles, 
including the photon, are spinless), where it forbids the 
pair creation diagram. Indeed, since the plus-component 
of the pair momentum is always positive, the pair can not 
be created by a virtual photon with q + = 0. As a result, 
when q+ — ► 0, the corresponding phase space volume 
tends to zero, and the amplitude of the pair creation dia- 
gram disappears. For systems involving fcrmions and/or 
vector photons, the amplitude of the pair creation dia- 
gram becomes indefinite if ui-q exactly equals zero, since 
it is given by an integral with an infinitely large inte- 
grand and zero phase space volume. For this reason, one 
has to take uj-q ^ 0. We set u>-q = a(u>-p), where a is a 
constant which we take positive, for definiteness. We will 
see below that for the rotationally non-invariant cutoffs 
discussed in Sec. UH a non-zero contribution to the EMV 
from the pair creation diagram survives, and moreover, 
it tends to infinity when a — » 0. 



ond order of perturbation theory, < Hf^Hemi >, is 
also of order eg 2 , it does not result in irreducible di- 
agrams. The Hamiltonian (|42p produces therefore the 
five contributions to the EMV, shown in Figs. [2HH] The 
triangle and pair creation diagrams are generated by 
< Hf bl H em i >, the left and right contact terms come 
from < HfbiH em 2 >, while < H em z > is responsible 
for the double contact term. Applying the rules of the 
LFD graph technique [l[ to these diagrams 1 , we can find 
analytical expressions for the corresponding amplitudes. 
However, before writing them down, one should note the 
following. When calculating the form factors in covariant 
LFD, the condition uj-q = on the momentum transfer 
is usually imposed (equivalent to q+ = qo + q z = in the 




FIG. 5: Double contact term 



1 The rules incorrectly prescribe to use the theta-function 8(u>-k) 
for the contact term. This theta-function should be removed. 



We can now proceed to the calculation of the LF di- 
agram amplitudes. The contribution of the triangle dia- 
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gram on Fig. [2] reads 

n 2 

p(trt) _ y I a, , i. \ a, iM . / 



(2tt) 3 



9{to-k)S(k 2 - ^)cf fc 



x (/ - # + ^t' + m)0[u;V - fc)] 
x<5[(p'-fc + cjT') 2 -m 2 ]^- 

x lp{i> - ft + ybr + m)9[u>-(p - k)] 
dr 

x 5{{p -k + lot) 2 - m 2 ] — (44) 



= WrJ dk± i & 

(/ — ft + 0t' + m)j p (p' — ft + i/jt + m) 
2u-(p' -k)r' 2uj-(p-k)T ' 

We have introduced here new integration LF variables 
and x in a standard fashion (see Sec. Ill A ip . The 
singularities of the integrand at x — and x = 1 are 
excluded by introducing an infinitesimal positive cutoff 
e. The values of r's are found from the conservation laws 
imposed by the delta- functions. Namely, 

2tu-(p — k)r = to 2 — (p — k) 2 

= (l- X )(^±^ + ^±jl-mA, (45a) 
V 1 — x x J 



2u-(p' - k)r' = m 2 - (p' - k f 



(l-x') 



m 2 k /2 L + /i 2 



1 -x' 



m 2 , (45b) 



where x' = oj-k/uj-p' — xj (1 + a), k' ± = — x'A, and A 
is the part of the three- vector q, transversal to us. Note 
that 



Q 2 (l + a) - a 2 m 2 



(46) 



The formulas (|45|) — (|46|) follow from the kinematical rela- 
tions listed in Appendix [A"l 

The contribution of the pair creation diagram, Fig. [3J 
is given by 



p(pazr) 
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(uj-k)S(k 2 - {i z )(Pk 



p (2tt) 3 

x (/ - ft + l6t' + m)9[uj-(p' - k)} 

dr' 

X6[(p' - k + ojt') 2 - to 2 ] — 

T 

x 7p(^ - ft + </>t' - tfon + m)6[u)-(k - p)} 
x5[(p - k + ujT 1 - ljti) 2 - to 2 ] — . 

Tl 

The line carrying the four-momentum ki corresponds 
to an antifermion and is described by the propagator 
(to — ft 1 )9(uj-ki)S(k 2 — to 2 ). It is convenient to introduce 



a new variable r 

„2 



-p(pai 
P 



: t — Ti , instead of Ti : 
f d 3 k 

dr' 

x0[u-(p' - k)]S[(p' -k + ujt') 2 - to 2 ] — 

T 



ft + i/it' + to) 



ft+ </rr + m)6[w'(k -p)] 
dr 



><7pI 

xS[(p - k + ujt) 2 - to 2 



(r' - r) 



(47) 
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d 2 ki 



l+a-e 



dx 



(2tt) 3 J ^J 1+t 2x 
(p" - ft + ybr' + m)-i P {i> - ft + i/it + to) 
X 2w(p'-fc)r' 2w(p-fc)(r'-r) ' 

The "shifted" r enters the delta-function and is de- 
termined by the same formula (|45ap as r in other 
T's. We can therefore use the same kinematics. The 
integration over dr by means of the delta-function 
6[{p — k + ujt) 2 — to 2 ] brings the factor \uj-(p — k)\ in the 
denominator. We used that \uj-(p — k)\ = —uj-(p — k) in- 
side the integration interval 1 + e < a; < 1 + a — e. 

The contribution of the left contact term to the EMV 
is shown in Fig. IDJa). The corresponding amplitude has 
the form 

9 1 "(uj-k)S(k 2 - ^ 2 )d*k 



r (/ct) = 
p 



(2tt) 3 

x(/ - ft + l6t' + m)9[uj-{p - k)} 

dr' 

X5[(p' - k + LOT') 2 - TO 2 ] — 

r 



; 7p 



2w(p - fc) 



d 2 fc. 



(48) 



(2tt) 3 

1+Q ~ £ (tf'-ft+ </>t' + m)j p </> 



2x 2u)-{p' - k)r' 2uj-{p - k) ' 
For the right contact term, Fig.fJJb), we get 



p(rct) 

p 



(2tt) 3 



?(uj-k)6(k 2 - Li 2 )d 4 k 



2w{p' - k) 
x9[v(p - k)]6[(p 
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d 2 k y 



(2tt) 3 

1_e dx (Z>7 P 



lp{i> - ft + i/>t + m) 
- k + lot) 2 - m 2 }— (49) 



+ (ZlT + TO) 



2x 2u)-{p' - k) 2u-(p-k)r' 
Finally, the double contact term, Fig. [5l yields 

,2 



p(2ct) 

p " (2tt) 3 



d 2 fci 



2x 2u\p' - k) 7p 2uj-{p - k) ' 







(50) 
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The full EMV is the sum of all the five contribu- 
tions igU), l[37] ) -([5n i) . Note that the limits of integra- 
tions over dx are different in these contributions. In or- 
der to make the calculations easier, we split the whole 
region of possible values of x into the three sub-regions: 
(1) e < x < 1 - e; (2) 1 + e < X < 1 + a - e; (3) 
l+a+e<x< +oo. We then represent each of the ver- 
tices Tp ot ' and T p 2ct ^ as a sum of integrals over these inte- 
gration sub-regions, removing the singularities by means 
of the same cutoff e: 



After some rearrangement of the contributions, we can 
represent the full EMV as 



(51) 



-p(Ici) _ p(l, let) , p(2,Zct) 
*■ p p ' P ' 



with 



where 



p(l,Zct) 

p 



p(2,ict) 

p 



1 



d 2 k± 



64tt 3 (wp) 2 
" 1_e dx — ft + </>t' + m)"f p i/) 
x 

2 1 



X) T 



64tt 3 (wp) 2 
,>1+Q - C 

X 



1 + a- x)(l 
d 2 k ± 

• - ft + i/jt' + m)^ p (/} 



(1 



x){l 



r(A) 
p 



p(c) 
p 



p(tri) , p(l,ict) , p(rct) . p(l,2ct) 
p P P P 



p(B) _ r (pair) , p(2,ici) _ p(l) 



p(2,2ct) , p(0) , p(l) 
L p P ^ P ' 



The two new integrals, 



r(o) 

p ' 



(52a) 
(52b) 
(52c) 



and, analogously, 



where 



p(l,2ct) 

p 



p(2,2ct) 
P 



p(2ct) _ p(l,2ct) 
P 1 P 



ff 2 ^7p ^ 
64tt 3 (lu-p) 2 

r-l— e 

X 
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+ 



p(2,2ct) 
P 



fix 



x(l + a — x)(l — x) ' 

^7p f J2, 



64tt 3 (wjj) 2 

x / 

^7p ^ 



crfci 



fix 







64tt 3 {wpf 



x(l 
fi 2 /fc± 



x)(l 



dx 



x(l 



x)(l 



p(o) 
p 



-(1) _ 



^7p^ 

647r 3 (l + a) {uj-p) 2 

g 2 tfiiptfi 



647r 3 a {uj-p) 2 



d z k. 



d 2 k. 



dx 

x 



1+a-e 



dx 



(53a) 



x(l + a — x) 

(53b) 



have been introduced in order to make the integrands 



of T p A ^ and T)^' non-singular in x, although each term 
on the right-hand sides of Eqs. (|52a[> and (|52b[) contains 
logarithmic divergencies either at x = or at x = 1, or 
at x = 1 + a. 



Taking the sums (|52ap and (|52b[) . we easily find 



P 



P 



16tt 3 



d 2 k. 



647r 3 (w-p) 2 



dx 
o x 



d 2 k A 



ft + m)"f p (p' — ft + m) 



^7p(^ 



(^ - 2p'-fc)0 2 ~ 2 P- fc ) 4(1 + a)(w-pf 



1+a 



dx 



x(l + a — x) 



- ft + </>t' + m)j p (p' —ft+ t/br' + m) <fil P fi 
(x — 1)t'(t' — t) a 



(54a) 
(54b) 



Since the regularization in x is no more required, we set The remaining part of the EMV, T p , is given by a 

6 = 0. 
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sum of regularized integrals 
9 2 filptfi 



p(C) = 
p 64tt 3 {uj-pf 



d 2 k A 



1 



l-e 



1 

+ 
■ 



1+a+e 



l+e 



1 + a 
dx 
x(l — x) 

dx 



dx 
x 



x(l + a — x)(l — x) 



(55) 



5. Electromagnetic form factors 

We represent the EMV via the following decomposi- 
tion: 



u'T p u = v! |fi7 p 



zF 2 
2m 



(p + y) p -2 7p 



B 2 - 



uj-p 



B. 



m 2 (Jjujp 
' (u-p) 2 



u, 



(56) 



which is similar, although not identical, to the one used 
in Ref. [2(J. We shall come back to this difference in 
Sec. IIVI The decomposition is determined by the five 2 
form factors Fi 2 , -Bi-3- We shall call F 12 the phys- 
ical form factors, in contrast to the non-physical ones, 
Bi—3, which must be absent in the physically observed 
EMV. The appearance of the three extra form factors is 
a property of LFD. Note that they appear in standard 
non-covariant LFD as well, but in this approach they can 
not be separated out from the physical form factors, and 
an illusion may occur that the EMV structure is deter- 
mined, according to Eq. (j3"fj)) , by the two form factors, 
as in the Feynman case. This is however not so, because 
the electromagnetic current operator in LFD has five in- 
dependent matrix elements, not two. One may argue 
that even if each particular LF diagram produces a con- 
tribution of the form to the EMV, the extra three 
structures disappear (i.e., one has -B1-3 = 0) after sum- 
ming up all the contributions, at least in a given order 
of perturbation theory. It would be so if the amplitudes 
of the LF diagrams were given by convergent integrals, 
which happens for the pure scalar case (scalar bosons 
plus scalar "photon" and no fermions). For systems in- 
volving fermions, however, the amplitudes strongly di- 
verge. Their regularization, as will be shown below, may 
give rise to the appearance of extra (w-dependent) spin 
structures in the perturbative and non-perturbative reg- 
ularized EMV's. In this sense the situation is analogous 
to the case of the fermion self-energy discussed in the 
previous section. 



2 The coincidence of the number of form factors with the number 
of LF diagrams is, of course, by chance. 



The procedure to calculate the form factors is quite 
similar to that exposed in Ref. [2(| • We define the matrix 



(p" + m)T p {p' + m) 



Am 2 

and the following contractions: 



(57) 



ci 

C2 
C3 

c 4 
C5 



Tr{O p7 p }, 

iq v 

2m 

(p + p'Y 



Tr {O p 0}, 



uj-p 



(58) 



LU-p 



■Tr{Op}, 
'°;Tt{O p9 b}. 



Now, using Eq. (|56|) . we can get the following relations 
between the form factors and the quantities Ci_s: 



ci 

C2 

C3 

C4 
C5 



(2 - 4r?)Fi - 677F2 + 2(477 + a)Bi 

+ (2 + a)B 2 + 2{l + a)B 3 , 

677F1 + 2r7(2 - ?7)F 2 - 277(4 - a)Bi 



-r?(2 + a)B 2 



-B 3 , 



2(2 + a)F 1 -2r 1 (2 + a)F2 
+4[2t7(1 + a) + a]Bi + (2 + a) 2 B 2 
+2(2 + o)(l + a)B 3 , 
(2 + a) [Fi - r/F 2 + aB x ] , 
a 2 

2(1 + a)F x - —F 2 + 2a(l + a)B u 



(59) 



where 77 = Q 2 / (Am 2 ). Solving the system of linear equa- 
tions (|59"j) with respect to Fi :2 and Bi_3, we express them 
through ci_5. These expressions are lengthy and we do 
not give them here. Then, substituting the formulas (|54p . 



p 



and r p C ' into Eq. ([5T|) . calculating 
the traces (I55|) , and expressing the scalar products of the 
four- vectors p, p' , k, and uj through the two-dimensional 
vectors kj_, A, and the scalars x, a, we cast each form 
factor in the form of a three-dimensional integral over 
d 2 k±dx. The expressions for the scalar products through 
the integration variables and the momentum transfer are 
given in Appendix El 



4- Regularization with rotationally non-invariant cutoffs 

We give here the final expressions for the form factors 
found by using the rotationally non-invariant cutoffs 
and e imposed on the variables |kj_| and x. As we said 



above, the vertex functions T p A ^ and T^' do not require 
to introduce a cutoff in x. As far as V p is concerned, 
it is represented by a sum of integrals, each requiring 
regularization and thus depending on e. However, as can 
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be established by direct integration in Eq. 
itself has a finite limit when e — > 0: 

r (C) = 5 2 Al log(l + a) (Z>w p 



55l). the sum 



32tt 2 a(l + a) (w-p) 2 ' 



(60) 



We therefore remain with the cutoff A± only. Note that 
the mutual cancellation of the terms singular in x, which 
happens for the full EMV, is connected with rather weak 
(logarithmic) divergence of the corresponding integrals, 
so that the cutoff Aj^ is enough to make them finite. 
The same took place in the fermion self-energy case (see 
Sec. Ill A2|) . As we shall see below [Eq. (|C11|) in Ap- 
pendix [C] , the quadratically divergent term ~ fx>u> p k\ 

results also from T p B ^ (but not from ) and cancels in 

(c) 

the sum with T p . 

The details of calculations of all the contributions 
(A,B,C) are given in Appendix [CJ Adding all them 
together, we arrive at the following final result for the 
form factors: 



Fa 

F 2 
Bi 
B 2 
B 3 



A i 



1ft 2 l0 g- 

lD7T z m 
2n2 



9 

4tt 2 



log- 



9 Z Q Z A m _ 9 
247r 2 m 2 V ° g M 8 



i 7 
7 ~ 8 

0(Q 4 ), 



g 2 Q 2 



3g 2 

167r 2 327r 2 m 2 

„2 



0(Q 4 ), 



64tt 2 ' 

„2 



327T 2 ' 

g 2 (m 2 



\Qir 2 m 2 



-log- 



A 



9 



32n 2 m 2 



— +A»M 2 log- 
a \ fi 



(61a) 
(61b) 
(61c) 
(61d) 

. (61e) 



Note that the terms ~ log ^ cancel in the form factors 
B\ and B 2 . These form factors are finite and do not 
depend on any cutoff, in contrast to B3 which does not 
have a finite limit when Aj_ — >■ 00. As we have already 
mentioned, neither of the form factors depend on the 
cutoff e. At the same time, the amplitude of each of the 
LF diagrams shown in Figs. (SHU diverges like A^ log e. It 
means that the senior divergent terms ~ A^ log e and ~ 
A^ cancel after the incorporation of all the LF diagrams. 
However, the form factors i?i_3 which must be absent 
in the physical EMV remain non-zero values. The same 
happened for the term in the full self-energy, Eq. 0. 

We also see that B 3 has a pole at a = 0. One cannot 
set here Q 2 = 0, since it is impossible to keep fixed a — 
(u)-q) / (u)-p) when q — > 0. More precisely, since A 2 > 0, 



from Eq. (|46|) we have Q > . However, nothing 
prevents us to take a — > at fixed Q 2 . As the inspection 
shows, the singular term ~ 1/a in B3 results from Y 1 ^" 11 ^ , 
Eq. g7D. 

The standard renormalization procedure (I37|) affects 



the form factor F\ only, so that 
9 2 Q 2 



1 



log- 



(1 



0(Q 4 



(62) 



The cutoff-dependent term ~ log(Aj_/m) cancels, as it 
should. The renormalized form factor Kp™ coincides 
with F2, Eq. (|61bp . It is important to note that in or- 
der to eliminate the three extra form factors -B1-3, one 
should introduce new w-dependent counterterms into the 
interaction Hamiltonian. Since B3 non-trivially depends 
on q through Q 2 = —q 2 and a = uj-q/ui-p, its decomposi- 
tion in powers of q contains an infinite number of terms. 
Therefore, to eliminate B3 we need an infinite number of 
local counterterms, for this particular type of regulariza- 
tion. 



5. Invariant Pauli-Villars regularization 

We can now calculate the form factors within LFD, 
but using the invariant PV regularization instead of the 
LF cutoffs A^ and e. 

Regularization by a PV boson. We start with the reg- 
ularization by one PV boson only. Since each of the LF 
diagrams contains one bosonic propagator, this regular- 
ization reduces to a simple subtraction 



PV, b 



T{Q 2 , m, (j,) - f(Q 2 , m, 



As before, T denotes any of the form factors Fi or Bi 
In the limit of large /Ji we have 



Fi(Q 2 ,m, /Lti) 
F 2 (Q 2 ,m,fj, 1 ) 



16tt 2 
0, 



log 



A, 



16tt 2 



, Mi 1 

l0g T 

m 4 



for arbitrary finite Q 2 . Subtracting these expressions 
from those given by Eqs. (I61a|) and (I61bj) . respectively, 
we get 



PV, b 



9 



167T- 



log 



/'l 



9 

At: 2 



log- 



9 2 Q 2 



log- 



F. 



PV, b 



247T 2 TO^ \ [I 

3g 2 g 2 Q 2 



16tt 2 327r 2 m 2 



n 15 
~ ~ 16 

f 0(Q 4 ), (63a) 
0(Q 4 ). (63b) 



In spite of the fact that Eq. (|63a|) differs from Eq. (|61a| , 
both equations lead to the same renormalized form factor 
F£ en , Eq. (|6"2"1) . The renormalized physical form factors 
obtained for the invariant and non-invariant types of reg- 
ularization do therefore coincide. 

Concerning the non-physical form factors, the situa- 
tion is quite different. Since Bi_3 were calculated for 
arbitrary /i, we find from Eqs. (|61c|) - (l61c| : 



B 



PV, b 



B. 



PV,b 



0. 



(64a) 
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b pv, b = 9 2 {p\ ~ M 2 ) j A± 
167r 2 m 2 (l + a) m 

.... ., ., i Mi - M 2 - V log - - 2 M f log ^ 

o27T 777 \ fj, m 

(64b) 

The formulas (|63|) are valid for /j — > and /ii — > oo, while 
Eqs. (|64[) hold for arbitrary and fix- The non-physical 



form factors B\ 2 turned into zero, while B PVb ^ and 
still diverges logarithmically. 

Regularization by one PV boson plus one PV fermion. 
Since the PV regularization by one boson only does 
not cancel the form factor B3, we have to introduce in 
addition a PV fermion. In that case, no contact terms 
appear at all, and the only diagrams which contribute 
to the EMV are the triangle (Fig. [2]) and pair creation 
(Fig. [3]) diagrams. Each of them contains two fermion 
propagators, both being subject to the PV regulariza- 
tion. For this reason we need to know the vertex with 
different internal fermions. We denote the masses of 
these fermions by 777^ and m^, with the indices i and i' 
being either or 1. Let 777,0 = rn an d 7771 be the physical 
and PV fermion masses, respectively. Analogously, we 
denote the physical and PV boson masses by fio = A 4 
and fix. 

The PV regularized EMV is defined by 



v pv, b+f 
p 



E (- 

i j> ,j—0 



iy+t +Jr p (777i,777j/,/J i ), 



and fi — > fij. The vertex T p , Eq. (|55[) . turns into zero 
already by the boson PV regularization and does not con- 
tribute to the PV regularized form factors. 



It can be shown that we do not need to introduce the 
PV fermion to regularize the form factors Fx. 2 and B\2- 
In other words, in the limit 777.1 — > 00 we would obtain the 
same formulas (f63|) - (|64al) as without the PV fermion at 
all. So, the only thing to do is to calculate B PV ' b+ ^ . This 
calculation is in principle quite similar to that performed 
above, but the algebra is more lengthy. We represent 
B PVb+f as follows: 



B 



where 



PV,b+f _ 



= E E B 



(N)PV,b 



{mi,mv), 



i,i'=0 



N=A,B 



(65) 



Bi N)pv > b (m 2 , 777,0 = Y^i-lYB^imumv,^) (66) 



3=0 



where the notation T p (mi,mi> , fij) stands for the cor- 
responding initial EMV calculated for the internal par- 
ticles with the masses m„ 777^, and fij. The quantity 
Tp(mi, , fij) itself is given by a sum of Eqs. (154 
changing everywhere 



77! 
77! 



p f 1 '"t ! 
777 2 — (p — k) 2 

2uo-{p-k) 1 
m 2 , — (p' — k) 2 
2u)-{p'-k) : 



and B% (to, , rrii' , fij ) is found through the vertex 

Tp(mi,mi',fij). Though the order of summations in 
Eqs. ([65|) and (|66[) does not matter, it is convenient to 

calculate first B^ N ' FV ' 1 '(m,-,r?v), separately for the two 
vertices, then add the results, and finally sum over the 
PV fermion indices. Omitting rather tiresome manip- 
ulations, we give here the expressions for the functions 
, rrii' ) , in order to see the details of the can- 
cellations: 



where 



B. 



(A)PV, b 



{mu rrii') = - 



g 2 (n 2 -fi 2 ){2 + a)^A ± 



167T 2 777 2 (1 + a) 2 
,2 



log 



9 



777 

/'1 . ■ ■ rn 



2f4 log ^ + 2fi 2 log 3{fi(- fi z 

rn fi 



B 



[B)PV,b 



327T 2 777 2 

,2/, ,2 ,,2\ A „2/,,2 ,,2 



R{rrii,mi'), 



{rrii, m^) 



g 2 {ni-fi 2 ) Ax , ST (pi ~ M ) m , 
log 1 — - — - R{mi,mi>), 



167r 2 777 2 (l + a) 2 



7/i 



167r 2 777 



(67a) 
(67b) 



R{rrii, rriii) 



64n 2 m 2 Q 2 



o o , rU7 „n 777,777,/ 

rn 2 - rn 2 log -f - 2Q 2 log 

mi 



SQ log 



(m 2 -m 2 ) 2 -(Q 2 -s Q ) 2 
(m 2 - rn 2 ,) 2 - [Q 2 + s Q ) 2 
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and s Q = y/mf - 2to 2 (to 2 - Q 2 ) + (m 2 + Q 2 ) 2 . Both 
contributions (|67|) are rather complicated functions de- 
pending on the physical and PV masses, as well as on 
Q 2 . In their sum, however, the dependence on the 
fcrmion mass m\ and on Q 2 drops out completely, and 
f° r ^2n=a b B^ N ^ PV ' b {mi, m,/) we arrive at the same ex- 
pression (|64b[) found previously without any fermion PV 
rcgularization. The final summation over the fermion PV 
indices, as prescribed by Eq. (|65[) . turns -B3 into zero: 



B 



PV, b+f 



0. 



(68) 



We can thus formulate the main results of this section: 
the full cancellation of the three non-physical form fac- 
tors is achieved by the double (i.e., bosonic + fermionic) 
PV regularization, whereas the single boson PV regular- 
ization cancels only two of the three non-physical form 
factors. The renormalizcd physical form factors obtained 
by using the PV regularization (both single and double) 
coincide with those calculated through the rotationally 
non-invariant cutoffs. 



B. Calculation in the four-dimensional 
Feynman approach 

The amplitude of the Feynman triangle diagram which 
determines the EMV in the given order of perturbation 
theory is 



r 

p (2tt) 4 



d 4 k 



1 



[k 2 -fi 2 + i0] 
'■ — ft + m)j p ($ — ft + m) 



[(p - k) 2 - m 2 + i0] [{p' - k) 2 - m 2 + iO] ' 



(69) 



where p 2 — p' 2 — m 2 . As we shall see, the spin structure 
of the Feynman EMV depends on how it is regularized, 
in contrast to the self-energy, Eq. (fT9|) . In order to avoid 
overloading the formulas by additional indices indicating 
the type of the cutoffs, we use in this section the same 
notations for the Feynman EMV and the form factors as 
for the LFD ones. 



we get 



*9 / ,2 



32tt 4 



d k±dk + dk 



i' — # + m)j p 



[{p' + - k+)(p'_ - &_) - (p' ± - k ± ) 2 -m 2 + iO] 
+ to) 

X [(p+ - k+)(p- - fc_) - (pi - k ± ) 2 - to 2 + iO] 

x [k+k-~ki- v 2 + ioy (70) 

We take the reference frame where p = 0. We then 
choose p + = p' + . The latter condition is equivalent, in 
the language of covariant LFD, to a = (u-q) / (w-p) = 0. 
In our LFD calculation of the form factors in Sec. IIII A[ 
we initially kept a to be non-zero and went over to the 
limit a — * after summing up all the LF contributions. 
A smooth limit a — > was important in the LFD frame- 
work, since the amplitudes of different diagrams shown 
in Figs. [2HS1 depend on a. The Feynman amplitude (|6"9"]) 
is well-defined for any values of its arguments p and p' , 
in particular, for p + = p' + . We can therefore safely set 
p + — p' + = from the very beginning. 

Similarly to Eq. |22|) for the self-energy, we represent 
the full EMV as a sum of the two terms: 



r, 



p(p+a) , -p(zm) 



(71) 



corresponding to the " normal" contribution and the zero 
modes, respectively. 

Like the self-energy case, the "normal" part T p p+a ^ is 
determined by the sum of the pole and arc contributions. 
It is calculated in Appendix iDl Sec. ID 21 The result reads 



p(p+a) 
P 



9 



16tt j 



d 2 k. 



1 dx 



(72) 



(/i 2 - 2p'-k)(^ 2 - 2p-k) A{uj-pf 



Comparing the right-hand sides of Eq. (|72p for Yp P+a ^ 

and Eq. (|Ma|) for the part T p A) of the EMV calculated 
within LFD, we see that they exactly coincide with each 
other at a = 0: 



Regularization with rotationally non-invariant cutoffs 



As we did for the self-energy, Sec. Ill B 1\ we first 
find the form factors by integrating the Feynman am- 
plitude (f6"9"l) in terms of the LF variables restricted by 
the transverse and longitudinal cutoffs. We shall see that 
such a procedure generates extra (non-physical) form fac- 
tors £?i_3 [see Eq. ([5^1) ]. in a very similar way as the use 
of the LFD rules. 

Rewriting the integrand in Eq. (|69j) through the plus-, 
minus-, and transverse components of the four-vector k, 



-p(p+a) 
P 



HA) 
p 



a=0 



(73) 



We have already encountered a similar situation above, 
in Sec. Ill B II where it was shown that the pole plus 
arc contribution to the Feynman expression for the self- 
energy exactly coincided with the full LFD self-energy 
[see Eq. (l24l) ]. Evidently, the vertex T p p+a ^ can be also 
represented via form factors, in the form of the decom- 
position ([56]) . but with its "own" form factors F^ a \ 



B 



(p+a) 
1-3 



Form factors for Tp^ were found in Appendix [Cl 



Due to the identity (|T3"|) . in order to derive the form fac- 



tors F^~ a> and B^^ 11 ' , we can make use of Eq. (|C2[) with 



,(p+a) 
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the coefficients from Eqs. (IC7|) - (IC9|) . taken at a = 0. We 
thus find 



? (p+a) 



9 2 , &± 

lD7T z m 



g 

4tt 2 



log- 



9 2 Q 2 f, m 



F, 



Cp+o) 



3.9 2 



O+o) _ 



16tt 2 327r 2 m 2 

~ 2 . A ± 
log 



.9 



9 

~ 8 
0(Q 4 ), 

„2 



0(Q 4 



(74a) 
(74b) 
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(p+o) 



B 



O+a) _ 



16tt 2 



327r 2 m 2 
5 2 



64tt 2 



m 
(6r, 



32tt 2 



(74c) 
MQ 2 ), (74d) 



327r 2 m 2 
-(4to 2 - 2^t 



4/x 2 + Q 



21og — + 

•Q 2 MQ 2 )] , 



log 

m 



(74e) 



where the function <p(Q 2 ) is given by Eq. IjClOp . Note 
that ip(0) = 0. As before, the form factors are 
calculated for \i — > and decomposed in powers of Q 2 , 
whereas -B^_1^ are exact in this sense. 

The zero-mode contribution Pp 21 ™-* is calculated in Ap- 
pendix |DJ Sec. ID 31 The expression for r p zm ^ is given by 
Eq. (|D2ip . We represent it in the form of decomposi- 
tion [|56p. Then the form factors in this decomposition 
have the form 



if 

b[ z 



(zm 



(zm 
2 

.(zm 



92 log^ + TrA^ 2 ), 



16tt 

9 2 , A ± 
log — 



r/v 



16tt 2 



8tt 



8tt 2 



y(Q 2 ), 



' /J (6m 2 -4 A( 2 + Q 2 )log^ 
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(75a) 
(75b) 

(75c) 

(75d) 
9 2 



32-7r 2 m 2 



32tt 2 
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A* ( 2 log - 



1 



(4m 2 - 2^ 2 + Q 2 MQ 2 ) 



(75e) 



and 175 



Ft 

F 2 
Bi 
B 2 
B, 



2n2 
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5 /, ™ 



.9 2 Q 



247r 2 m 2 
3ff 2 



log- 
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,2n2 



16tt 2 327r 2 m 2 
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0(Q 4 



64tt 2 ' 

„2 



32tt 2 ' 

„2 



32tt 2 



(76a) 

(76b) 
(76c) 
(76d) 
(76e) 



The expressions for £?i_3 are exact, i.e., these formulas 
are valid for any Q 2 and /i. That is, the Q 2 -dependence 

coming from B[ p _+ a \ Eqs. ([73c ]H [75e ]) . as well as the Ax- 
dependent terms are exactly canceled by the correspond- 
ing zero mode contributions b[ z ^ . The form factors 
Fi t 2 and B\ 2 are the same as those obtained from the 
LFD diagrammatic approach, in Eqs. (|6T|l . while B^s in 
Eqs. (|61ep and (|76ep . are different. Moreover £?i_3 are 
not zeroes, though we started from the Feynman expres- 
sion for the EMV, which initially had no any relation to 
the light front (i.e., to the four- vector u>). We will discuss 
this situation in more detail below, in Sec. IIVI 



2. Invariant regularization 

To complete our analysis, we give here the form fac- 
tors obtained from Eq. (|69p by using an invariant regu- 
larization of the divergent integral over cftk. The stan- 
dard recipe consists in using the Feynman parametriza- 
tion, then making the Wick rotation and imposing a cut- 
off A 2 on the modulus of the Euclidean four-momentum 
squared, similarly to how it has been done in Sec. Ill B 21 
After that, identifying the structure of the EMV fS9"]) 
with the decomposition (|36p . we reproduce the well- 
known result: 3 



F 2 = 



9 



32tt 



A 2 g 2 / to 15 

log — + log — 

m z 47r z \ 11 16 



9 



24ir 2 m 2 
3.9 2 



log- 
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~ 8 
0(Q 4 ). 



0(Q 4 ), (77a) 



(77b) 



167T 2 327T 2 TO 2 

Eq. (|77bp exactly coincides with Eq. (|76bp obtained 
for the rotationally non-invariant cutoffs. The for- 



The final expressions for the form factors are given by 
the sum of the corresponding quantities from Eqs. I|74p 



3 ^MQ 2 = 0) is an anomalous magnetic moment. Introducing the 
coupling constant a = cj 2 /(47r), we get in the Yukawa model: 
F 2 (0) = 3a/(4?r) that differs by a factor of 3/2 from the QED 
value a/ (2tt). 
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mulas (|77a|) and (|76a|) differ by a constant (i.e., Q 2 - 
independent) part. The corresponding form factors F[ en 
however also coincide after the renormalization (|37| . 

In principle, one might identify the structure of the 
Feynman EMV (regularized by means of an invariant 
cutoff!) with the LF decomposition (|56|) containing five 
form factors. In this case the formulas (f77[) would be 
found for F± t 2, while one would arrive at the evident re- 
sult -B1-3 = for the other form factors. 

The regularization of the EMV (|69|) by means of one 
PV boson results in the same formulas ([77]) . changing A 2 
to 



IV. 



DISCUSSION 



We have shown above that the perturbative spin- 1/2 
fermion self-energy and the EMV calculated in LFD de- 
pend on the orientation of the LF plane, when the tradi- 
tional regularization in terms of the transverse and lon- 
gitudinal cutoffs is applied to the corresponding ampli- 
tudes. This dependence reveals itself in the appearance 
of extra spin structures in the decompositions of the self- 
energy and the EMV in invariant amplitudes. In covari- 
ant LFD on the plane ui-x = with w 2 = 0, it is con- 
veniently parameterized through the four- vector in. The 
corresponding decompositions for the self-energy and the 
EMV are given by Eqs. (JTJ) and (|55|) . respectively. The 
structure of the LF self-energy is characterized by the 
three scalar functions A, B, and C = C(p 2 ) + Cf c , in- 
stead of the two usual ones, while the on-energy-shell 
EMV contains five form factors, the two standard F\ 2 
and the three extra, Bis, ones. 

Performing calculations of the self-energy by integrat- 
ing the Feynman amplitude written in terms of the LF 
variables with the same cutoffs as in LFD, we did not 
encounter any m-dependence in the final result, Eq. (|29[) . 
The latter therefore does not coincide with the LFD self- 
energy (UJ. For the EMV, both methods lead to in- 
dependent structures with five form factors. However, 
the EMV calculation by means of the LFD rules results 
in the set of form factors (|6Tj) , whereas the calculation 
of the Feynman amplitude in terms of the LF variables 
gives a different set of form factors f76|) . Though this 
difference concerns the form factor B3 only, the corre- 
sponding EMV's do not coincide nevertheless with each 
other. 

The extra form factors do not, of course, appear when 
the Feynman amplitude is calculated with a spherically 
symmetric cutoff in four-dimensional space. This can be 
achieved, e.g., by imposing a direct cutoff on the Wick 
rotated integration variable in Euclidean space \k 2 \ = 
k 2 + k\ < A 2 or by using the PV regularization. In the 
case of LFD, we deal with three-dimensional integration 
variables. Therefore, constructing rotationally invari- 
ant (i.e., ^-independent) cutoffs which restrict a three- 
dimensional integration domain encounters serious diffi- 
culties since the integration domains in LFD amplitudes 



differ from each other, and it is not easy to restrict them 
simultaneously in a self-consistent way. This is a rea- 
son why the PV regularization looks much more prefer- 
able. It is naturally generalized to LFD and allows us 
to remove ultra-violet divergencies in an w-independent 
way. We showed that in order to cancel completely all in- 
dependent terms in the LF self-energy and the EMV it is 
necessary to use the PV subtractions for both the boson 
and fermion propagators. Note that this cancellation oc- 
curs for arbitrary (finite) PV masses. Simultaneously, all 
ultra-violet divergencies disappear. After that, the reg- 
ularized LFD results coincide exactly with the Feynman 
ones calculated with the same PV subtractions. 

The physical quantities like the functions A and B in 
the self-energy or the form factors F\ z in the EMV, can 
be easily extracted from the corresponding amplitudes. 
Once we know the self-energy £(p) explicitly, taking the 
traces 



Tr{S(p)} - Ag 2 A(p 2 



Tr{S(p)(n} = 



Ag 2 {uj-p) 



B(p 2 ) 



allows to get A{p 2 ) and B(p 2 ), because the term ~ i/b in 
the self-energy ([Jj) does not contribute to these traces. To 
separate the physical electromagnetic form factors F12 
from the non-physical ones, Bis, it is enough to con- 
sider the plus-component of the current (or, in terms of 
covariant LFD, the contraction of the current with tn p ) 
written in the form (|56p . Indeed, in the limit ui-p — m-p', 



u'T + u = u (in-r)u = u 



iF 2 
Ira 



u. 



As we see, Bis dropped out from here. This is how- 
ever not always the case since the decomposition (l56l) is 
not unique. Take, for instance, the decomposition of the 
same EMV given in Ref. [2CJ |: 



u T p u = u < Fi7 p + —a pv q v 



1 



_l_ 

oj-p m(l + rj) 



(p + p')p 



LU-p 



B. 



m 2 (/]U) p 
(uj-p) 2 



u, (78) 



where 77 = Q 2 /4m 2 . It is easy to see that the form factors 
£?i_3 in Eqs. (|56p and ([75)1 are identical, while 



F =Fi 



2r?gi 
I + 77' 



2Bi 
l + rj' 



If we now identify Fi t % with the physical form factors, 
they will differ from F1.2, unless B\ = 0. If we used for 
the regularization the LF cutoffs A^ and e, we would get 
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Bi = -g 2 /64w 2 [see Eqs. ([BTc| and (|76c)) ] and 



integral turns into 



- g 2 , A_l g 2 ( m 7 
F l = T7TJ> lo S + 739 I lo S 



16tt 



fa 



47T 2 \ ' fJ. 

m 21 

24-7r 2 m 2 y ° fjb 16 
7 ff 2 5 ff 2 Q 2 



log- 



8 

0(Q 4 



32tt 2 1287r 2 m 2 



o(g 4 ). 



do not coincide neither with the Feynman form fac- 
tors, given by Eqs. (|77p (if we identify with A), nor 
with the LFD ones defined by Eqs. (|76aj) and ([76b|) . The 
regularization by the PV boson kills B\ and this ambigu- 
ity disappears. However, if the non-physical form factors 
are not canceled completely, the result is sensitive to the 
form of the EMV representation. 

Although one can use the plus-component of the cur- 
rent ([56]) to extract the physical form factors in hadron 
phenomenology, it is not enough in many other cases, 
when the knowledge of the full matrix structure of the 
vertex is needed. For example, this occurs in LFD non- 
perturbative approaches, when the EMV enters as an off- 
energy-shell subgraph into a more complicated diagrams. 
In such a case, one can not simply ignore the extra spin 
structures, because they may contribute to observable 
quantities. The latter statement concerns the self-energy 
as well. 

As our results show, the dependence on the LF plane 
orientation appears not only in LFD amplitudes, but also 
in Feynman ones, if we put the cutoffs in the LF vari- 
ables and x, i.e., constrain the integration domain 
by a spherically non-symmetric region, the orientation of 
which follows the orientation of the LF plane. At first 
glance, this seems contradictory: (i) the extra form fac- 
tors -B1-3, [Eqs. (|76c|) — (|76e|) ] originate from the spheri- 
cally non-symmetric cutoffs, but (ii) they do not depend 
on the cutoff values. This is in fact a normal situation 
which can be illustrated by the following toy example. 
Consider the three-dimensional integral 



kj d k 



(k 2 +' 



2)5/2 



Fib 



(79) 



It is logarithmically divergent at infinity: F\ ~ 
Jd|k|/|k|. Let us introduce the spherically symmetric 
cutoff |k| < L. A simple calculation gives 



Fl -| log ^T^ + T (log8 - 4) + (z 



Let us now regularize this integral by a cutoff imposed 
on the two-dimensional variable k^ in the plane or- 
thogonal to an arbitrary direction n. That is, we set 
k^ = k 2 — (k-n) 2 < Aj_. In other words, we integrate 
over the volume of a cylinder of the radius Aj_ and of in- 
finite length, with the axis directed along n. The initial 



jreg 



k<i kj 



(k 2 



[Ai - k 2 + (k-n) 2 ] d 3 k 
(80) 



Fidij + Bi(Sij - n^). 



The integral along n (or over dk z , if n is parallel to z) 
converges, similarly to the convergence of the integral 
over dx in the Feynman amplitude written through the 
LF variables. We see that the spherically non-symmetric 
cutoff Aj_ generates one extra form factor B±, like this 
happens for the Feynman amplitude. We obtain 



Ft 
Bi 



2?r 



log- 



Ai 



O 



27T ^ ( 1 

T + AT 



The leading terms oc logi in F\ and F\ coincide, pro- 
vided we identify Aj_ with L. However, there is a dif- 
ference in the finite parts: i*i — F± = 4^ (log 8 — 4), like 
the difference between the term j| in Eq. (|77a[) and | 
in Eq. (|76ap . The transverse cutoff A_l generates a finite 
extra form factor B\, which itself does not depend on 
Aj^, similarly to -Bi_3 in the EMV. This example clearly 
mimics the properties of the above calculation of the elec- 
tromagnetic form factors from the Feynman amplitude 
regularized by the LF cutoffs. 

It is therefore misleading to think that one can de- 
rive some "true" LF amplitude starting with the covari- 
ant Feynman amplitude and calculating it in the LF 
variables: / = J...d k = j . . . d 2 k±dxdk-. This in- 
tegral diverges (except for some particular cases) and 
it has no sense without regularization. It depends on 
the size, shape and orientation of the integration domain 
constrained by the regularization procedure. In other 
words, there is no "covariant Feynman amplitude" in it- 
self. There exists only an inseparable couple: the co- 
variant Feynman amplitude together with the rules to 
regularize it. 

When we calculate amplitudes by means of the LFD 
graph technique rules, the non-invariant integration do- 
main is not the only source of the w-dependence. An- 
other source is hidden in the rules themselves. Indeed, 
these rules are not obliged to reproduce exactly the am- 
plitude which follows from the Feynman approach, even 
for the same integration domain D. The divergent (and, 
after regularization, cutoff-dependent) terms are treated 
non-identically in the two approaches and may be differ- 
ent. This fact is illustrated by the above calculation of 
the form factor B$ (cf. Eq. (|61e[) obtained within LFD 
with Eq. (|76ep coming from the Feynman approach) . The 
renormalization procedures may be different too. How- 
ever, the renormalized, observed amplitudes must be the 
same. 

When the w-dependent contributions to the self-energy 
and EMV survive, one has to cancel them with appro- 
priate extra counterterms in the interaction Hamiltonian. 
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These counterterms are inherent to the LFD Hamiltonian 
and should not be confused with the traditional charge 
and mass counterterms in the original Lagrangian. The 
need for such a counterterm for the renormalization of 
the self-energy was already advocated in Ref. Q, where 
the non-perturbative fermion mass renormalization was 
studied in the two-body approximation within covariant 
LFD. This new counterterm cancels the contribution ~ 1/1 
in the self-energy ((7}. 

However, the introduction of additional specific coun- 
terterms may help only if their number is finite and 
tractable. For the case of the fermion self-energy we 
would need only one w-dependent counterterm, since the 
coefficient C is a constant. In the case of LFD calcu- 
lations of the EMV, the form factor B3 depends non- 
trivially on the momenta (through Q 2 and a), which 
would generate an infinite number of local counterterms 
to kill it, even in perturbation theory. For this reason, the 
use of counterterms can not serve as an universal tool for 
the calculation of rcnormalizcd quantities in LFD with 
the traditional transverse and longitudinal cutoffs. 

If the renormalization is done correctly, the depen- 
dence of any renormalized amplitude on the cutoffs dis- 
appears and the final result is the same, regardless to the 
type of the regularization used. However, from the prac- 
tical point of view, the difficulty of calculations rapidly 
increases with introducing extra counterterms. There- 
fore, rotationally invariant regularization (like the PV 
one) seems by far more preferable in the LFD framework. 
The approach developed in Ref. [21| . regularizing field 
theory by defining fields as distributions, could provide 
another method of rotationally invariant regularization. 



nian, which also depend on the LF plane orientation. 
Taking them into account, one can calculate the renor- 
malized amplitudes. 

The use of spherically symmetric (in four-dimensional 
space) regularization, like, for instance, the PV one, con- 
siderably simplifies calculations. In perturbation theory 
it allows to avoid the presence of these extra countert- 
erms at all. 
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APPENDIX A: KINEMATIC AL RELATIONS 

We give below some kinematical relations used in the 
self-energy and EMV calculations within covariant LFD. 
Amplitudes of the LF diagrams are expressed through 
the three-dimensional integrals over d?k/Ek. Take the 
four- vectors k = (ek,k) and p = (pcbP)j and construct a 
new four-vector R = k — xp with x = (uj-k)/ (uj-p). Since 
ijj-R = and oj 2 = 0, in the reference frame, where p = 0, 
we have £k — xpo — |k|| | = 0. Taking into account that 

J 



= Jkl + k| + /i 2 and changing p to 
at the relations 



we arrive 



V. CONCLUSION 

We have calculated perturbatively the fermion self- 
energy and the EMV in the Yukawa model, in two differ- 
ent ways: (z) by the LFD graph technique rules, taking 
into account all necessary diagrams; (ii) by integrating 
the Feynman amplitudes in terms of the LF variables 
(i.e., the transverse and longitudinal parts of momenta) 
and summing up the pole, arc, and zero-mode contribu- 
tions. For the same set of the cutoffs imposed on the 
LF variables, both methods give different results for the 
regularized amplitudes. This difference disappears if the 
invariant PV regularization is used instead of the rota- 
tionally non-invariant one. 

Such properties follow from the fact that the cutoffs 
constrain a spherically non-symmetrical integration do- 
main, the symmetry being destroyed by the choice of a 
distinguished direction defined by the orientation of the 
LF surface. The dependence of regularized amplitudes 
on the LF plane orientation is conveniently taken into 
account in the explicitly covariant version of LFD by con- 
structing extra spin structures. To exclude these struc- 
tures from the physically observed quantities one may 
introduce new counterterms in the interaction Hamilto- 



Iky I = 



2 2 
x p 



2x\/p 2 



£k = 



k\ + (i 2 + x 2 p 2 
2x\fp 2 



In the variables kj_ and x, the invariant phase space ele- 
ment becomes 



d 3 k/sit; = d 2 k±dx/x. 



(Al) 



We can now express the scalar product k-p entering the 
quantity r, Eq. (j4|), through k^ and x. For this purpose, 
we represent the invariant quantity R 2 in two different 
ways: 



and 



R z 



R 2 



(k - xp) 2 = 1? - 2x(k-p) + x 2 p 2 



(e k - xp ) 2 
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Equating the right-hand sides of these expressions yields 



k-p = 



x 2 p 2 



2x 



(A2) 



In the calculation of the EMV we have one more four- 
vector p' . We define two new four- vectors R = k — xp and 
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R 1 = k — x'p', where x' = (uj-k) / '(tv-p') = x/(l + a). From 
the equalities lu-R = uj-R' = 0, in the reference frame, 



where p = 0, follows R 



-(k 
Eq. 



Ri 



R 



/2 



R 



/2 



— x'A.) 2 , where A = = q^. 
~2\), we get from here 



Analogously to 



k-p = 



k-p' 



+ p? + x 2 m 2 
2x ' 
(kj. - x'A) 2 + fi 2 + x' 2 m 2 
2x~> 



(A3a) 
(A3b) 



From these formulas follow the expressions (I45aj) 
and (|45b[) . For completeness, we give here the other 
scalar products which are needed for the EMV calcu- 
lations: 



oj-k 


— x(uj 


■P), 




ui-p' 


= (n 


- ol)(uj 


■P), 


p-p' 


_ Q 2 
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+ m 2 




CO 2 


- o, 


P 2 = 


P = 



, K = \i . 

After integrating over the azimuthal angle of the vector 
k^, the form factors depend on A 2 , and one should relate 
it with the invariant square of the momentum transfer 
Q 2 = —{p 1 — p) 2 . If one had co-p = uj-p' (equivalent to 
a = 0), it would be simply Q 2 = A 2 . In our case the 
relation is more complicated. Indeed, since oj-p — u>om, 
we have, on the one hand, uj-p' = (1 + a)uJom, and, on 

the other hand, uj-p' = w (y / A 2 + p'jj + m 2 - We 
thus find 



IP II 



- m 2 a(2 + a) 
2m(l + a) 



and 



Q 



Pi 



-my 



A 2 + a 2 m 2 



From here Eq. (|46[) follows. 



(A4) 



to get the results in a form closer to LFD, we will use the 
residues. The integrand has two poles in the points 



and 



fc_ = k^ = p- 



k 2 ± + (i 2 - iO 



xp + 



iO 



(1 - x)p^ 



If — oo < x < — e or 1 + e < x < +oo, both poles lay, 
respectively, above or below the real axis. If e < x < 1 — e 
the poles are situated on the opposite sides from the real 
axis. We close the integration contour by an arc of a 
circle of the radius L, either in the upper half-plane or in 

(2) 

the lower one, so that the pole k_ ' is always outside the 
contour. The situation is illustrated in Fig. [51 where the 
integration contour is shown for the case e < x < 1 — e. 

Note that the order of integrations we have chosen in 
Eq. (|20|) requires some care in treating the cutoffs. Since 
the integration over dk- is performed first, we should 
keep L arbitrary large, while the other cutoffs e and A^ 
are considered as being finite. In other words, we imply 
that L ^> A 2 L /(me). Such a convention ensures mutual 
position of the poles and the contour, as exposed above. 
Then £ p+a (p) is represented as 



km 




-L I 


/ +L 




J 




arc. 



FIG. 6: Contour of the integration over dk- 
poles are shown for e < x < 1 — e. 



Positions of the 



APPENDIX B: DERIVING THE LIGHT-FRONT 
SELF-ENERGY FROM THE FEYNMAN 
APPROACH 

Proceeding from the Feynman amplitude (|19p for the 
self-energy, Sec. IIIBI we calculate here the pole, arc, and 
zero-mode contributions. 



1. Pole and arc contributions 

We first consider the integral over dk- in Eq. (j2"0)) for 
the x-integration in the limits — oo < x < — e, e < x < 
1 — e, and 1 + e < x < +oo, determining Y. p+a (p). It can 
be calculated directly (by primitive). However, in order 



jg 2 p^ 

32tt 4 



d 2 k\ 



^arc, low dx 



1+e 



, (Bl) 



where S po ; e is the residue at the pole k^\ multiplied by 



-27ri, while £ 



arc, low 



and E. 



arc, up 



come from the integra- 
tions along the arcs in the lower and upper half-planes, 
respectively. The result for E po ; e reads 



-'pole 



2-7ri(^ — fi + m) 



p+\k 2 L + m 2 x — p 2 x(l — x) + n 2 (l 



(B2) 



The four-vector k here is an on-mass-shell four-vector 
(k 2 = pi 2 ) with the components expressed through k^ 
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and x as follows: k = (fe_, k+, kj_) = ( 
To calculate S arc / oto and X„ 



xp A 



^arc,up, we should move along 
the arc in the clockwise and counter-clockwise directions, 
respectively. In the points of the arc fc_ = Le 1 ^ and 
dk- = iLe z ^d(f>, where is an azimuthal angle. In the 
limit L — > oo, we retain in the integrand the dominating 
/c_-term. We thus obtain 



Taking the sum of (IB5|) and (|B6|) , we find that it has the 
form of eq. ([Ml): ^zm{p) = 9 2 C zm ^ with 



*2*'J rfkil ° g k 



2 ' \i 2 



Calculating this integral, we find Eq. (|27|) for C z 



7+ 



-'arc. Low 



2p\ x(l 



x) 



(B3) 



Substituting Eqs. jB2| and (JH3J) into Eq. (JHU) and going 
over to the explicitly covariant notations by means of the 
identities 7+ = yb, p+ = uj-p, we obtain that Yi p+a (p) is 
given by Eq. ([23]) . 



APPENDIX C: LIGHT-FRONT 
CONTRIBUTIONS TO THE FORM FACTORS 

In this Appendix, starting with the LF vertices T p A \ 
T ( P B \ and T P C) found in Sec. lIIIA21 Eqs. (|52aj) - (j52c|) . we 
calculate their contributions to form factors Fi 2, -Bi-3- 



2. Zero modes 

Consider now the zero-mode term H zm {p). We denote 
E zm (p) = sW(p) + EW(p), 

where the two items on the right-hand side correspond 
to the contributions to Eq. (|2H)) from the infinitesimal 
integration regions — e < x < e and 1 — e<a:<l + 
e, respectively. Take first T,{°m(p) (bosonic zero modes). 
If x = 0, the term k-p + x in the denominator of the 
integrand disappears and the integral over dk- diverges 
linearly. This infinite contribution should therefore be 
proportional to 8(x). We take x — > and keep in the 
numerator and in the denominator of the integrand in 
Eq. (|20p the leading fc_-terms (where k- is not multiplied 
by x) only. That is 



E<°> = 



ig 2 p+ 
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■h+k- 



[k-p + x - - fi 2 + iO] (— k~p+) ' 
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The double integral over dxdk- is of the type of Eq. (|E1| 
with b = (kj_ + [i 2 )/p + and h = (see Appendix |E} 
Using Eq. (|E2|) . we obtain, in the covariant notations: 



(w-p)i 1 
l0 § ,-2 7 „2 " - - - 



(B5) 



Changing the variable a; — > 1 — x, for £im(p) (fermionic 
zero modes) we similarly get 



roW 327r 3 (w-p) 
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1. Preliminary transformations 



Each contribution to a form factor T [T = Fi or Bj) 



from the vertices T 
an integral of the form 



(A) r (B) 



or rl can be written as 



T = 



dx 



dk\ 
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#/(ki,x,</>,Q 2 ,a), 



^ ( C1 ) 
where is the angle between the two-dimensional vec- 
tors k^ and A, / is a given function (/'s are different 
for different form factors), and the modulus of A is ex- 
pressed through Q 2 by Eq. ppj) . The limits x m i n and 
Xmax depend on which vertex is considered. 

The analysis of the dependence of the form factors on 
the cutoff Aj_ (when it tends to infinity) allows to repre- 
sent each of them as 



T = aiA^ + a 2 log — + a reg , 
m 



(C2) 



where the coefficients a±, 02, and gular func- 

tions depending on the particle masses, Q 2 , and a, but 
independent of Aj_. If a 7^ 0, no other divergencies 
(i.e., cutoff-dependent terms) excepting those listed in 
Eq. (|C2[) appear. Since the coefficients a\ and 02 deter- 
mine the dependence of the form factors on the cutoff, 
we will calculate them at arbitrary a, in order to see 
in detail how the cutoff disappears (if so!) after taking 
into account all the contributions T p A \ T p B \ and T p C ' to 
the full EMV and its subsequent renormalization. As far 
as the coefficients concerned, we will find them 

in the limit a — > 0, retaining non-vanishing terms only. 
This simplifies calculations a lot, but leads to the same 
final result as if a was non-zero. 

Although direct integrations over dip and (ik^ can al- 
ways be done analytically for arbitrary Q 2 and a, the 
results turn out to be very cumbersome, so that the final 
integration over dx (if it can be done analytically) takes 
too much time. To make the computations more effec- 
tive, we will proceed in the following way. We take the 
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initial expression for a given form factor in the form of 
Eq. (|C1|) and study the behavior of the integrand as a 
function of in the asymptotic region |kj_| — > oo. For 
this purpose we decompose it in a Laurent series: 



/(kj_,a;,0,Q 2 ,a;) = f {x,Q 2 ,a) 

/i(x,Q 2 ,a) coscj) f 2 (x,<f), Q 2 ,a) 
|kx| + k 2 ± 



(C3) 



We use this decomposition in order to define the func- 
tions /o and f 2 (the term with f\ drops out after the 
integration over d(f>). We then represent / as a sum 



/(kj_,a;,0, Q 2 ,a) 



/re S (ki,x,0, Q 2 ,a) 
fl n 2 \ i <f>2(x,Q 2 ,a) 
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where 



1 f 27T 

$ 2 (x,Q 2 ,a) = — / f 2 (x,<t>,Q 



,a)d(j) (C5) 



and /3(ie) is a positive function which can be chosen in 
any convenient way in order to avoid the singularity at 
— > 0. Eq. (|C4[) should be considered as a definition 
of the function /^(k^, x, cf>, Q 2 , a). The latter, after the 
integration over dcf>, decreases faster than l/k^ in the 
asymptotic region, and its integration over dk^ does not 
require any regularization. The functions fo(x, Q 2 , a) 
and $2(2;, Q 2 , a) are rather simple, so that the two last 
addenda on the right-hand side of Eq. (|C4j) can be easily 
integrated over all the variables. 

After separating out the terms which slowly decrease 
when |kj_| — > 00, we remain with a set of regular func- 
tions freg- At arbitrary a, these functions are even more 
complicated than the initial integrands /. However, we 
calculate the regular contributions to the form factors in 
the limit a — > 0. So, it is enough to consider the limit- 
ing functions / reg (k 2 L , x, 4>, Q 2 , a — > 0) which are much 
simpler than /(kj_, x, <f>, Q 2 , a). If f reg is singular at 
a — ► 0, we decompose it in powers of a and retain all 
non- vanishing terms. 



The trick exposed above allows to find the coefficients 
in Eq. (JC2J. Indeed, substituting Eq. jC4| into Eq. (|CT|) . 
we get 
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dx 



dki 



x I d(j) f reg (k L ,x,(j), Q ,a -> 0) 

- 2ir<S> 2 (x, Q 2 ,a -» 0) log/?(x) J . (C6c) 

The integrations over and dk^ can be done analyti- 
cally for all the five form factors. For F\ f2 the results of 
the remaining integrations over dx at arbitrary Q 2 and 
/i are not expressed through elementary functions. In 
order to simplify the formulas, we decompose both F\ 2 
in powers of Q 2 up to terms of order Q 2 and take the 
limit /j, — ► 0. Concerning the form factors -Bi_3, all the 
integrations are performed analytically without any ap- 
proximation. 

Below, in Sees. IC 2"MC 41 we list the coefficients aj., a 2 , 
and a reg entering Eq. ()C2|) . pointing out in superscripts 
which form factor the given coefficient belongs to. The 
coefficients are calculated by means of Eqs. (|C6|) . sepa- 



rately for the vertices T p 
and (1551)1. 



(A) r (B) 



and r^ C) [see Eqs. 



2. Contribution to the form factors from Tp" 4 ' 

It is convenient to set in Eq. (|C4[) /3(x) = x 2 + ^(l— x). 
The limits of the x-integration are x m i n = 0, x max — 1. 
Then 
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(C8) 



16tt 2 (1 + a) 
„2 



.9 



167T 2 TO 2 (1 + a) 



(2 + a)ii 2 



{l + a)(2 + a)z a 
(3 + 3a + a 2 )m 2 Q 2 [(l + a)Q 2 - a 2 m 2 



1 + a 
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where z a = (2 + 2a + a 2 )Q 2 + 2a 2 m 2 , and 

.2 / m 7 



9 



Fi 

7 

reg 

reg 

i B2 

reg 
reg 



log 



o 2 Q 2 /, m 9 
log 



res 4tt 2 V M 8/ 247r 2 m 2 \ /i 8 
3 3 2 ff 2 Q 2 



o(g 4 ), 



16tt 2 327r 2 m 2 

-1 + 4^(Q 2 )] 



0(Q 4 



64tt 2 



(C9) 



32tt 2 

~2 



9 



32ir 2 m 2 



[i + ¥0 2 )] 



Ai 2 (21og^ + l) - (4m 2 - 2 M 2 + Q 2 MQ 2 ^ 



r 



with 

<p(Q 2 



I 4m 2 

V 



VQ 2 + 4m 2 + VO 5 



2m 



3. Contribution to the form factors from F 



(B) 



(CIO) get 
I 



We set = 1, x 



1, Xtj 



rain — x > ^max 



1 + a. We thus 



a : = a x = a t ^ = 0, 
g 2 log(l + a) 



327r 2 m 2 a(l + a) ' 



16tt 2 



1 - 



2(l + a)Q 2 + a 2 (2-a 2 )m 2 



a 9 = — 



(1 + a) 2 z Q 
5 2 (2 + a) [(1 + a)Q 2 + a 2 m 2 
16tt 2 (1 + a) 2 z a 



-F 2 _ 



16tt 2 (1 + a) 

„2 



a 9 = — - 



4a 2 (2 + 2a + a 2 )m 2 + [8 + a(2 + a) (8 + 2a + a 2 )]Q 2 
(l + a)(2 + a)z a 
2 (2 + a)m 2 Q 2 [(l +a)Q 2 -a 2 ™ 2 ] 



167r 2 m 2 (l + a) : 



1 + a 



r 



(Cll) 



While integrating over dx in Eq. (|C6cj) . we introduce 
a new variable £ = (x — l)/a. Then 

l+a pi 

dx(...) = a / #(...)• 



After that the integrands / reg should bedecomposed in sition ((SBJ). For this reason T p > contributes to B 3 only. 



Contribution to the form factors from F 



(C) 

The matrix structure of the vertex T p is the same as 
the one in front of the form factor B3 in the decompo- 
se) 



powers of a up to terms of order 1/a. Performing this 
transformation, we obtain the following result: 



From Eq. l[rJD|) we easily find 



"reg 

l Bl = 
reg 



<lg = 0, 



-a B2 

2 reg 



16-7T 



MQ 2 ), 



(C12) 



g 2 log(l + a) 



reg 327r 2 m 2 
with <p(Q 2 ) given by Eq. (fCTtlf 



^ + (4m 2 -2 M 2 + Q 2 MQ 2 ) 
a 



327r 2 m 2 a(l + a) 



all the other coefficients being zero. 



(C13) 
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APPENDIX D: DERIVING THE LIGHT-FRONT 
ELECTROMAGNETIC VERTEX FROM THE 
FEYNMAN APPROACH 



Proceeding from the Feynman amplitude ([69)1 for the 
EMV, Sec. IIII Bl we calculate here the pole, arc, and 
zero-mode contributions. 



1. Preliminary transformations 

Under the condition p' + — p + = 0, the LF components 
of the four- vectors p and p' are 



e<x<l — e, l + e<x< +oo and two infinitesimal 
ones — e < x < e, 1 — e < a; < 1 + e. The integrals over 
the "normal" regions are calculated by summing up the 
residue and arc contributions, as in Sec. Ill B II The con- 
tributions (if any) from the two infinitesimal regions of 
the integration over dx, which do not vanish in the limit 
e — > 0, correspond to the zero modes. Because of rather 
weak (logarithmic) divergence of the initial integral ([69)1 
we expect that the result of the full integration over dk_ 
and dx has a finite limit at L — > oo and e — ■+ 0, while 
the remaining integration over d 2 k± produces logarith- 
mic dependence of the EMV on the cutoff Aj_. 



P_L =0, p- = , 

Pi = q±, v- = — ; 

p+ 



(Dl) 



where q_j_ is the transversal part of the three-dimensional 
momentum transfer: 

qx=p±-px, Q 2 = -(p'-p) 2 = qi- 

Introducing a new variable x by means of the relation 
k + = xp + and using the notations 



ki+^ 2 



k^+m 2 , (kj_-q ± ) 2 +m 2 



P+ P+ 
we cast Eq. (P70)) in the form 

„2 



P+ 



(D2) 



r 



p 128ttV 



tfkt 



dx 



+ 

+°° dk_ [k- 1+ + M'h P [k- 1+ + M] 
viv 2 v' 2 



where 














Vl = 


xk- 


- a + i0. 










V2 = 


(x- 


l)k- - b 


-(x- l)p- - 


FtO, 






v' 2 = 


(x- 


l)fc_ - b' 


- {x - l)p'_ 


+ *0, 




(D4) 


M = 


(x- 


1)p+7- ~ 


P-1+ - 2k ± 


•7_l - 


2m, 


(D5) 


M' = 


(x- 


l)P+7- - 


P-1+ ~ 2(k 


l q 


l)"7± 


— 2m. 



The matrices Ai and M! do not depend on fc_. 

The integral in Eq. (|D3|) diverges at |kj_| — > oo. If 
x = or x = 1, it diverges also at fc_ — > ±oo. We con- 
strain the integration over dk 2 j_ in d 2 k± = ^dk^dcf) by 
the cutoff Aj_, and introduce also the cutoffs L in fc_, as 
in Eq. (|2ip . and e in x, splitting the region of the integra- 
tion over dx into three "normal" regions — oo < x < — e, 



2. Pole and arc contributions 



We start with the calculation of rp P+a - ) . The procedure 
is quite similar to that exposed in Appendix IB II and in 
the papers 0, E3, EE IlE El ■ The integrand in Eq. (|D3|) 
has three poles: 



fc (i) = a-iO 



(2) _ b+ (x - l)p- - iO 
ar-1 ' 

(3) _ b' + {x- l)p'_ - iO 
x-1 



k 

kf> = 



(D6) 



If — oo <i<-eorl + £<i< +oo, all of them are, 
respectively, either in the upper half-plane of fc_ or in 



the lower one. If e < x < 1 — e, the pole k_ is in the 

(2) (3) 

lower half-plane, while the poles k_ and k_ are in the 
upper one. We can thus write, analogously to Eq. (|B1[) : 



d 2 k 



128^ 

r arc,low dx 



(r p p ole -r a p rc ' low )dx 



T arc,up dx 



, (D7) 



where T p ° le equals the residue of the integrand in the 
pole fc_ = k^\ multiplied by —2-kv. 

%irip 2 + (p" - ft + m)7 p (^ - jt + m) 



-ppole 
P 



x [{p 1 — k) 2 — m 2 ][(p — k) 2 — m 2 ] 



(D8) 



with k 2 = [i 2 , and the arc contributions are 



F 



arc, low 



pare, up 
P 



■7+7p7+- (D9) 



x(l — x) 2 

Substituting Eqs. (JD8]) and (0 into Eq. ([DT]) . we find 
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p(p+a) 
P 



16tt 3 



d 2 fc. 



dx (f? — ft + m)j p (fi — ft + m) 
x \{j>' — k) 2 — m 2 ] [{p — k) 2 — m 2 ] 

I 



7+7 P 7+ , 1 
4p± e 



(D10) 



We substitute here log(l/e) by J dx/x and take the limit 
e — > 0. The latter is achieved simply by setting e = in 
the integration limits, because the integrand is no more 
singular neither at x = nor at x = 1. We also return 
to the covariant notations, replacing 7+ by 1/1 and p+ by 
w-p. In this way we find that r p p+a ' ) is given by Eq. (|72p. 



In order to calculate these integrals, we need their asymp- 
totical limit when L — > 00 and e — > 0. As explained in 
Sec. Ill B 11 one should take the limit L — > 00 first, while 
keeping e finite, and then allow e — > 0. 

We represent ^p X_0 ^ as 



3. Zero modes 

Let us now consider the zero mode contribution to 
Eq. (ID3p . It can be written as 



p(zm) 
P 



where 

£(*=0) = 



^p 



'.<T 



128tt 4 ^ 



d 2 k ± 



(Dll) 



L dfe_ [fe_7+ + >l , ]7p[fc-7+ + M] 

L ViV 2 v' 2 

(D12a) 

L [A:_7+ + X']7 p [/s_7+ + M] 



j(x=0) 



dx / <^ 7+7p7+ 7 



fc_ 1 
+ (A4'7 P 7+ + 7+7p.M) T +M'j p M T 

ViV 2 V 2 ViV 2 V 2 



For further calculations, we will make use of the formulas 



k n DJx) EJx) FJx) 

nv i , jiw^ n = 0, 1, 2, (D13) 



Vl 



l'2 



v i v 2 v 2 where v%, v 2 , and v 2 are defined by Eqs. (|D4|1 . and the 

(D12b) functions D n (x), E n (x), and F n {x) are 



A>(z) = 777 > W = 

1)^) = -^, E 1 {x) = 
t\t x 



x - 1 



D 2 {x) = — , £ 2 (z) 



ht 2 ' 

b+ (a - l)p- 
ht 2 ' 
[b+(x- l)p_] 2 



F (x) = 
Fi{x) = 
F 2 (x) = 



x-l 



b' + (x~ iy_ 
t[t 2 

[V + (x - l)p'_] 2 

(x - l)t[t 2 



(D14) 



r 



with 

t\ = (x — l)(xp- — a) + bx, 

t[ = {x-l){xp'_-a)+b'x, 

h = b-b' + (p- -p'_)(x-l). 

The values of a, b, and b' are defined by Eqs. (ID2[) . 

After the transformation (|D13[) . the problem of find- 
ing Gp X ~°^ reduces to the calculation of integrals of the 



following three types: 



dxf(x) 
dxf(x) 
dxf(x) 



dk_ 
-l vi 
L dk_ 

-L V 2 

dk- 
-l v' 2 



(D15) 



with various functions f(x), being either D n (x), E n (x), 
and F n (x), or their products with M and M.' . It is easy 
to see that the latter two integrals in Eqs. (|D15[) always 
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give zero in the limit e — > 0, while the first one is of the 
type of the integral (|E4|) from Appendix [E] Using the 
formula (|E4|) and taking into account that D 2 (0) = 1, we 
get 



C(*=o) 
^p 



= 7+7p7+ 
= -7+7p7+ 



T , / s f dk- 
/ dxD 2 (x) / 

7- £ 7-l wi 



TT 2 + 2-Kl I log — — log - 

a e 



(D16) 



The calculation of Q p x ^ is similar, but the algebra is 
more lengthy. Changing the variable x — ^ x — 1 , we have 



: (x=l) 



dx I dk- \ 7+7p7+ ~—r 

' - L I uiu 2 u' 2 



(A^i7 P 7+ +7+7 P A^i) 



Jfc_ 



U\U 2 U 2 



+M' 1 7 p Mi 



1 



U\U 2 U 2 



where 



Ml = {x + l)fc_ — a + iO, 
U2 = xk— — b — xp- + iO, 



a,, 



xk- — b' — xp'_ + iO, 



(D17) 



Mi = xp+7_ -p_7+ - 2k ± -7j_ - 2m, (D18) 
Xi = xp + 7_ -p'_7+ - 2(k_L - q±)-7_L - 2m - 
We have now, instead of Eq. (|D13[) : 

k n _ D n (x + 1) , E n (x + l) F n (x + l) 



U\U 2 U 2 



Ul 



«2 



n = 0, 1, 2. (D19) 

Again, we encounter integrals of the types (1D15[) . v's be- 
ing changed by it's. The terms with u 2 and u' 2 in the 
denominators contribute only, while those with u\ dis- 
appear in the limit e —> 0, after the integration over dx. 
Using the formulas (|E4[) and (|E5|) from Appendix [El we 
obtain after some transformations: 

Q [ P X=1) = 1+1 pi + |tt 2 + 2th flog ~ - log i 

+27TiFl(k_L)} 

+ 2TTIJ3+ [X' lo7p 7 + + 7+7 P X 10] ^2 (k x ) , (D20) 

where 

■Mio = -P-1+ - 2k ± -7 ± - 2m, 

M' w = -p'_l+ - 2(k ± - q_L>7 ± - 2m, 

(a-b')(b-b')+bp'_-b'p- b 
+ (b - b'f g V' 

H ^ = Wibh) lo 4- 



In order to find the form factors ^ zm ) from the vertex 
T p zm \ it is convenient to cast the latter in an explicitly 
covariant form. For this purpose we introduce the four- 
vectors 



r 
r' 



(r_,r + ,r x ) = (p-,0, -kj_), 
(/ ,r' rl) = (p'_,0,q ± -kx). 



Then, taking the sum of Eqs. (|Dl6|) and ([D20|) , substitut- 
ing it into Eq. (|D11|) . and changing everywhere 7+ — » ^, 
p + — * cj-p, we arrive at the desired result: 



-p(zm) 
P 



9 



32tt 3 



(/ + H7p— + — 7p(/+ m) 
OJ-p oj-p 



H 2 Qt±)>. (D21) 



To find the form factors, one should substitute T y p 
into Eq. ([57[) . instead of r p , then find ci_5 by means 
of Eqs. (|58|) . and finally revert the system of linear equa- 
tions ([55]) taken for a = 0. We give here the result: 



(zm) 



= F, 



(zm) 



= 0, 



b 



(zm) 



— — — — f d 2 k 1 _b l (k 1 _), 



where 



6i(kj.)=lT a (k 
&a(kj_) = - 



6a(k±) = -2iT 2 (k ± ), 
[ffi(kj.) - 2(k ± -cLL + m 2 )JJ 2 (k±)] 



It is convenient to use the following representations: 



fli(kj 



zp_ + (1 - z)p'_ + a - 6' 
6z + 6'(l - z) 

a — 



ff 2 (k ± ) = — / 
P+ Jo 



az + 6(1 — z) 
dz 

bz + b'(l - z) ' 



which can be checked by direct integration. Calculating 
first the integrals over d 2 k± and then over dz, we find the 
zero mode contribution to the form factors, Eqs. ([75)) . 



APPENDIX E: CALCULATION OF THE 
ZERO-MODE INTEGRALS 



Finding the zero mode contribution to the self-energy 
and EMV requires calculation of the asymptotic value of 
the integral 



h = 



dx 



dk- 



xk- 



b — xh + iO 



(El) 
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at L — > oo and e — > (6 and /i are considered as finite 
quantities independent both of fc_ and x). As we shall 
see below, the real part of I\ is finite in this limit, so that 



Reii 



Re 
Im 



dx 



dk- 



xk- — b — xh + iO 

+ OC 

dk- 



27rlm 



lim / dx I 

r+oc 

X / rfy e i(xk--b-xh+iO)y 
J\ 

r+oo 

n / dy e - l{b - l0)v 

J A 

x J dxe~ lhxy 8(xy) 

n f +oc , sin6y 
= -2tt / dy 

Jo V 
= — 7r 2 sgn(6). 

The imaginary part of I\ is divergent in L and e, and we 
have to retain finite values of the cutoffs: 

Im/i = — 7T J dx j dk- S(xk- — b — xh). 

Using the identity 5(z) = S(—z), we can easily see that 
Im/i docs not change under the transformation b — > —b. 
It is thus legitimate to substitute b by \b\. Calculating 
the integral by means of the (^-function, we get 



(1*1 \ 
(-— dx r dx} 
J- e \x\ + Jj^^j 

A L 2 -h 2 nl \ 
-7T I log — — h 2 log e 1 . 



Im7i 



Collecting the results together and neglecting the terms 
of order (l/L) 2 and higher, we finally obtain 



h = -7r-sgn(6) - r. ( log -logi 



(E2) 



Let us now consider the integral needed for the calcu- 
lation of the zero-mode contribution to the EMV: 



/:r/ 



(E3) 



_ L xk- — b — xh + iO 
where the principal value prescription for the pole at 



x = is implied. The real part of I 2 is found by di- 
rect integration: 



Rc/ 2 = 



dx 



P.V 



dx 

log 



L xk- — b — xh 
xL — b — xh 



= 

eL 



xL + b + xh 

A N 3 

eL 



We have expanded the result in a series in powers of 
l/(eL), retaining the leading term only. Note that ac- 
cording to our prescription we should first tend L to infin- 
ity while keeping e finite, and then tend e to zero. Under 
such a convention, we obtain Rc/ 2 = 0. The imaginary 
part of I 2 reads 



Im7 2 



[ e dx f 

—ir / — / dkS(xk- — b — xh) 

J-e x J-L 

r dx r*L-»-*» 

J- e ^ J- V (?/) 

2irh 
~b~' 



-xL—b—xh 



Finally, 



h = 



2-rrih 



From the results obtained above it is easy to derive the 
following expressions for the generalized integrals 



h 



f dxf(x) J 



dk- 



L xk- — b — xh + iO 



(E4) 



7r^sgn(6) + 2-Ki ( log — — log - 



2-Kih 



fix) 



/(0) 



dk- 



/(0), 



xk- 



b — xh + iO 



(E5) 



7T sgn(6) + 2m (log — - log - 



L 



1 



/'(0), 



where f(x) is an arbitrary function supposed to be 
smooth and finite at x = 0. 
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